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ABSTRACT 



It has long been regarded as difficult if not impossible for a cosmological model to account 
simultaneously for the galaxy luminosity, mass, and velocity distributions. We revisit this is- 
sue using a modern compilation of observational data along with the best available large-scale 
O ' cosmological simulation of dark matter. We find that the standard cosmological model, used in 

*— > , conjunction with halo abundance matching (HAM) and simple dynamical corrections, fits - at 

^ ' least on average - all basic statistics of galaxies with circular velocities Vi^iic > 80 km s~^ cal- 

culated at a radius of ~10 kpc. Our primary observational constraint is the luminosity- velocity 
relation - which generalizes the TuUy-Fisher and Faber- Jackson relations in allowing all types of 
galaxies to be included, and provides a fundamental benchmark to be reproduced by any theory 
O^ ' of galaxy formation. We have compiled data for a variety of galaxies ranging from dwarf irreg- 

OO . ulars to giant ellipticals. The data present a clear monotonic luminosity-velocity relation from 

Cn ' ~50 km s~^ to ~500 km s~^, with a bend below '-^ 80 km s~^ and a systematic offset between 

late- and early-type galaxies. For comparison to theory, we employ our new ACDM "Bolshoi" 
VO ' simulation of dark matter, which has unprecedented mass and force resolution over a large cosmo- 

^P , logical volume, while using an up-to-date set of cosmological parameters. We use halo abundance 

matching to assign rank-ordered galaxy luminosities to the dark matter halos, a procedure that 
automatically fits the empirical luminosity function and provides a predicted luminosity- velocity 
relation that can be checked against observations. The adiabatic contraction of dark matter halos 
k^ , in response to the infall of the baryons is included as an optional model ingredient. The resulting 

}^ ' predictions for the luminosity- velocity relation are in excellent agreement with the available data 

5^ 1 on both early-type and late-type galaxies for the luminosity range from M^ = —14 to Af,. = —22. 

We also compare our predictions for the "cold" baryon mass (i.e., stars and cold gas) of galax- 
ies as a function of circular velocity with the available observations, again finding a very good 
agreement. The predicted circular velocity function is also in agreement with the galaxy velocity 
function from 80 to 400 km s~^, using the HIPASS survey for late-type and SDSS for early-type 
galaxies. However, in accord with other recent results, we find that the dark matter halos with 
T4irc < 80 km s~^ are much more abundant than observed galaxies with the same Vcirc- Finally, 
we find that the two-point correlation function of bright galaxies in our model matches very well 
the results from the final data release of the SDSS, especially when a small amount of scatter is 
included in the HAM prescription. 

Subject headings: cosmology: theory — dark matter — galaxies: halos — galaxies: structure 



1. Introduction 

The cosmological constant + cold dark matter 
(ACDM) model is the reigning paradigm of struc- 
ture formation in the universe. The presence of 
large amounts of dark mass in the surroundings of 
galaxies and within galaxy clusters has been es- 
tablished firmly using dynamical mass estimates 
that include spiral galaxy rotation curves, velocity 
dispersions of galaxies in clusters and x-ray emis- 
sion measurements of the hot gas in these systems, 
as well as strong and weak lensing of background 
galaxies. The ACDM model also correctly predicts 
the details of the temperature and polarization of 



the c osmic background radiation ([Komatsu et al. 
20101 ). A few issues remain where the model and 



the observations are either hard to recon cile or 
very difficult to compare ( Primackl 120091 ) . Ex- 
amples of this are the so-ca.lled r nissing satel- 



lites problem (iKlvpin etal.l 11999 : 



1999: Bullock et al.l 120001: iWillman et al 



Moore et al. 
2004t 



Maccio et al.l 120101 ) and the cusp/core nature 
of the central density profiles of dwarf galaxies 



(IFlores fc PrimacL 
19971: |Valenzuela_ 



In all three cases we need to estimate the circu- 
lar velocity (a metric of dynamical mass) at some 
distance from the center of each dark matter halo 
that hosts a visible galaxy. Unfortunately, theory 
cannot yet make accurate predictions for the cen- 
tral regions of galaxies because of uncertain bary- 
onic astrophysics. As a compromise, we propose 
to use the distance of 10 kpc. Measurements of 
rotational or circular velocities of galaxies either 
exist for this distance or can be approximated by 
extrapolations. At the same time, theoretical pre- 
dictions at 10 kpc are also simplified because they 
avoid the complications of the central regions of 
galaxies. 

Our LV relation is a close cousin of the Tully- 
Fisher (TF) relation and, indeed, we will use some 
observational results used to construct the tradi- 
tional TuUy-Fisher relation. However, there are 
substantial differences between the TF and the 
LV relations. The standard Tully-Fisher rela- 
tion tells us how quickly spiral disks rotate for 
given luminosity. The rotation velocity is typ- 



20101: Ide Blok 20101 ) 



a , " tu^ irrrmr; — Tt^^l — - ^^ ^ ncally measured at 2 .2 disk scale lengths (e.g., 
HiBllpllde B ok fc McGaugl fe ;^^^^;^^ ^^ ^^^^^ ^^^ ..^J, Jl^ 

a, et al.l 20071 : Cjovernato et al. /. , j u \ . -i . i x x- i 

~ (i.e., stars and cold gas) contribute a substantial 

fraction of the mass. Instead, at the 10 kpc radius 
used here for the LV relation, the dark matter is 
the dominant contribution to the mass in all but 
the largest galaxies. More importantly, the LV 
relation includes not only spiral galaxies, but all 
morphological types. Thus, the LV relation is the 
relation between the galaxy luminosity and the to- 
tal mass inside the 10 kpc fiducial radius. 

In order to make theoretical predictions for 
the LV relation, we need to estimate the lumi- 
nosity of a galaxy expected to be hosted by a 
dark matter halo (including those that are sub- 
structures of other halos). There are different 
ways to make those predictions. Cosmological iV- 
body-|-gasdynamics simulations will eventually be 
an ideal tool for this. However, simulations are 
still far from achieving the resolution and phys- 
ical understanding necessary to correctly model 
the small scale physics of galaxy formation and 
evolution. Early simulation s had problems repro- 
ducin g the TF relation (e.g.. iNavarro fc Steinmeta 
2000t ). lEke et al.l (|200l[) could reproduce the slope 



An outstanding challenge for the ACDM model 
that we address here is to reproduce the observed 
abundance of galaxies as a function of their over- 
all properties such as dynamical mass, luminosity, 
stellar mass, and morphology, both nearby and at 
higher redshifts. A successful cosmological model 
should produce agreement with various observed 
galaxy d ynamical scaling laws s uch as the Faber- 
Jackson (Faber fc Jacksonlll976r) and Tully-Fisher 
(JTuUv fc Fisheij|l977l) relations. 

Making theoretical predictions for properties of 
galaxies that can be tested against observations 
is difficult. While dissipationless simulations can 
provide remarkably accurate predictions of vari- 
ous properties of dark matter halos, they do not 
yet make secure predictions about what we actu- 
ally observe - the distribution and motions of stars 
and gas. We need to find a common ground where 
theoretical predictions can be confronted with ob- 
servations. In this paper we use three statistics 
to compare theory and observations: (a) the lu- 
minosity - circular velocity (LV) relation, (b) the 
baryonic Tully-Fisher relation (BTF), and (c) the 
circular velocity function (VF). 



of the TF relation, but created disks that were 
too faint by about 0.5 magnitudes in the /-band 
for any given circular velocity. Recently the situa- 



tion h as improved. For example, iGovernato et al 
(J2OO7I ) produced disk galaxies spanning a decade 
in mass that seem to fit both the /-band TF rela- 
tion and the baryonic TF relation very well, as 
well as the observed abund ance of Milky Way - 
type satellites. Most recentlv. lGuedes et al.l(|201ll ) 
have produced perhaps the best match yet to a 
Milky Way-type galaxy in ACDM using a high- 
resolution smooth particle hydrodynamics simula- 
tion with a high density threshold for star forma- 
tion. 

Making predictions for a large ensemble of sim- 
ulated galaxies is yet another challenge. Semi- 
analytical models (SAMs) are a way to make 
some progress in this direction. These models 
have the advantage of producing large-number 
galaxy statistics. They typically include many 
free parameters controlling the strength of the 
various processes that affect the build-up of the 
stellar population of a galaxy (i.e., cooling, star 
formation, feedback, starbursts, AGNs, etc.). Un- 
fortunately, these normalizing parameters can 
be difficult to c onstrain obscrv ationally (e.g., 
Somerville fc Primack 1999; Benson fc BoweJ2010 



The models aim to reproduce the observed num- 
ber distributions of galaxies as a function of ob- 
servables such as luminosity, stellar mass, cold 
gas mass, and half-light radius, along with scal- 
ing laws such as the TuUy-Fisher relation and the 
metallicity-luminosity relation. 

Early S AMs suffered from seriou s defects. The 
models of iKauffmann et al.l (|l993l ) were normal- 
ized using the observed TF relation zero-point, 
which resulted in a luminosity function with a 
very steep faint- end. On th e other hand, models 
such as those of ICole et al.l (|l994) were normal- 
ized to reproduce the observed "knee" in the LF 
but this resulted in a large offset in the TF rela- 
tion zero-point. Later models have shown mod- 
erate success in reproducing eit her the luminos- 



ity f unction (Benson et al.l 



tion ([Somerville fc Primacklll999l ). but it has been 



2003) or the TF rela- 



difhcult to match both simultaneousl y when ro- 



tatio n curves are treated re alistically (jCole et al 
2000D- iBenson et al.l (|2003l ) used a combination 



of disk and halo reheating to obtain reasonable 
agreement with the observed LF except at the 
faint end, where they still overpredict the number 
of d warf galaxies. If the WMAP 5-year cosmol- 
ogy (JKomatsu et al.ll2009l ) were used, their models 



would also produce too many very bright galaxies. 
The TF relation they obtain has the correct slope 
but their disks are too massive at any given lu- 
minosity. Most recently, iBenson fc Bowed (|2010l ) 
used a sophisticated version of their GALFORM 
semi-analytic model to obtain sets of parameters 
that minimize the deviations from twenty one ob- 
servational datasets including the LFs in several 
bands and at different redshifts, the TF relation, 
the average star formation rate as a function of 
redshift, clustering and metallicities among many 
others. Not surprisingly, even their best model 
has difficulty fitting such a large number of simul- 
taneous constraints. In particular, the LF in the 
K band overpredicts the number of dwarf galax- 
ies by almost an order of magnitude at the faint 
end, while the LFs at high redshift consistently 
overpredict the abundance of all galaxies. In ad- 
dition, the halos they obtain contain too many 
satellite galaxies, resulting in too strong a galaxy 
two-point correlation in the one-halo regime. The 
TuUy-Fisher relation of their best fit model also 
shows a systematic offset of about 20 — 40 km s~^ 
towards higher circular velocities for any given lu- 
minosity when compared to observations. 

Recent high-resolutio n A^-body co s molog i- 



cal simulations su ch as ISpringel et al.l (120051 ): 
Klvpin et al.l (J2010I ) have volumes large enough 
to obtain the mass function of dark matter 
(DM) halos, but there is no direct way to com- 
pare it to observational measurements of the 
luminosity or stellar mass functions of galax- 
ies. A new technique recently emerged that al- 
lows us to bridge the gap between dark mat- 
ter halos and galaxies. It is commonly re 



ferred to as abundance matching dKravtsov et al 



2004l:lTasitsiqmi et al.l20o3:IVale fc OstrikeiJ2004 : 



Conroy et al 



'2006!; 



Conrov fc Wechsleil 12009: 



Guo et al. 2010; Beh roozi et al.ll2010l ). Halo abun 
dance matching (HAM) resolves the issue of con- 
necting observed galaxies to simulated dark mat- 
ter (DM) halos by setting a one-to-one correspon- 
dence between red-band luminosity and dynamical 
mass: more luminous galaxies are assigned to more 
massive halos. By construction, it reproduces the 
observed luminosity function. It also reproduces 
the scale depende nce of galaxy cluste r ing over a 



range of epochs (jConrov et al.l l2006t iGuo et al 



2OIOI ). When abundance matc hing is used for th e 
observed stellar mass function ( Li fc Whitell2009l ) . 



it gives a reasonably good fit to the lensing results 
(JMandelbaum et al.1120061 ) on the relati on between 



the s tellar mass and the virial mass (IGuo et al 
2OIOI) 



Guo et al.l pOld ) also tried to reproduce 



the observed relation between the stellar mass and 
the circular velocity with partial success: there 
were deviations in both the shape and the ampli- 
tude. At circular velocities Vc = 100—150 km s~^ 
the predicted circular velocity was ^ 25% lower 
than the observed one. They argued that the dis- 
agreement is likely due to the fact that they did 
not include the effect of cold baryons. Below we 
show that this is indeed the case. 

The paper is structured in the following way. 
Section[2]describes in detail the observational sam- 
ples used to compare with the results of our anal- 
ysis. Secti on [3l briefly describes our new Bolshoi 
simulation (JKlvpin et al. 2010 ) and compares it to 
other large cosmological simulations. In section |4] 
we describe some characteristics of dark matter 
halos. Section [5] describes the abundance match- 
ing method used to relate observed galaxies to 
the DM halos in the Bolshoi simulation and ex- 
plains the procedure used to measure key quanti- 
ties such as the circular velocity for these model 
galaxies. Section [5] describes in detail the obser- 
vational samples used to compare with the results 
of our analysis. Section [6] shows the LV relation, 
the baryonic TuUy-Fisher relation, the galaxy cir- 
cular velocity function and the galaxy two-point 
correlation function obtained using our model and 
compares them to the observations described in 
Section [2j A brief comparison with related results 
in the literature is given in Section [T] Section |8] 
presents a discussion of our results and Section |9] 
summarizes them. 

2. Observational Data 
2.1. Late- type galaxies 

We use several datasets to construct the LV 
relatio n for observed galaxies 



Springob et al 



(J2007I ) compiled a template /-band TuUy-Fisher 
sample of 807 spiral galaxies of types Sa-Sd in or- 
der to calibrate distances to ~ 4000 galaxies in 
the local universe. Template galaxies were cho- 
sen to be members of nearby clusters in order to 
minimize distance errors. Their photometry con- 
tains distance uncertainties so the scatter should 
be taken cautiously and only as an upper limit to 



the intrinsic TF scatter. Circular velocities were 
obtained using HI line synthesis observations or 
optical Ha rotation curves when those were not 
available. The maximum circular velocity was ob- 
tained by using a model fit to the observed profiles. 
Since the authors correct for the effects of tur- 
bulence by subtracting 6.5 km s~^ linearly from 
the velocity widths, it was necessary to de-correct 
them by adding this term back in to obtain the 
true circular velocities. 
The 



Pizagno et al. 



([20071) sample was se- 
lected from the S loan Digital Sky Survey (SDSS) 
( York et al.ll2000l ). It is one of the most complete 
and unbiased samples available of Ha rotation 
curves of disk galaxies and was studied in an at- 
tempt to accurately measure the intrinsic scatter 
in the TF relation. Luminosities were taken from 
SDSS r-band photometry, yielding the best match 
with the luminosities assigned to our model galax- 
ies. For this sample we used the asymptotic value 
of the rotation velocity they obtained using a func- 
tional fit to the rotation curves. 

In order to test the predictions of the ACDM 
model using abundance matching (ACDM -t- HAM 
for short) with the largest dynamic range possible, 
we included in our comparison the l atest Tully- 



Fisher dwarf galaxy sample studied bv lCeha et al. 
(|2006l ). Their sample consists of about 110 late- 



type galaxies with luminosities measured in the 
r-band and rotation velocities measured using HI 
emission. 

The three samples above constitute our major 
observational dataset. We further cut them by 
selecting galaxies with high inclinations (i > 45° 
or axis ratio b/a > 0.7) to minimize uncertainties 
due to projection effects. Additionally, we include 
only galaxies with better than 10% accuracy in 
the measurement of the maximum circular veloc- 
ity. These cuts leave a total of 972 galaxies in the 
major sample. 

For comparison with the datasets mentioned 
above, we include other smaller samples found in 
the literature. The sample of Blanton ct al. (20Q8|) 
was also obtained from SDSS and is comprised of 
only isolated galaxies with hi gh inclinations. The 



HI galaxy sample used by Sa kai et al.l (|2000l ) 



was 

selected to have small scatter for use as a distance 
calibrator. It is important to note that while the 
fit shown here minimizes both the errors in rota- 
tion velocity and in luminosity, it may be artifi- 



cially shallow due to selection effects. 

Certain assumptions about galaxy colors had 
to be made in order to convert the different obser- 
vational samples to the "'^r-band measurements 
we chose for our model. In order to convert the 
J-ban d luminosities measured by ISpringob et al 



(|2007l ) to the r-band, we cross-referenced their 
data with the sample of IPizagno et al.l ( 2007 ) and 
used the median (r — /) colors of the galaxies 
present in both catalo gs. To convert fr om the 
i?-band magnitudes of ISakai et al.l ( 20001 ) to the 
SDSS "'^r-band w e used the tran sformation equa- 
tions obtained by iLuptonI (|2005t ) along with the 
typical °'^(r — z) co lor of disk ga laxies in the SDSS 
sample studied bv iBlanton et al. (2003a ). In ad- 
dition, f or redshift zero data se ts, the k-correction 
given in iBlanton et al.l ( 2003bl ) was used to con- 
vert from z = to z = 0.1 photometric bands. 

Lastly, since the obscuring effect of dust extinc- 
tion as a function of disk inclination is corrected 
for in TuUy-Fisher samples but not in observed 
LF estimates, we had to de-correct the luminosi- 
ties of the spiral galaxies in all of the TF samples. 
To do so, we estimated and added the median ex- 
tinction in the r-band as a function of rotation 
vel ocity using th e met hod and sample employed 
bv IPizagno et al.l (|2007i) . This correction is ^ 0.4 
mag for the brightest disks, declining to ^ 0.3 mag 
for 14irc « 100 km s~^. Th e se va lues are close to 
those found bv ITuIIv et al.l (|l998l ) for the extinc- 
tion in a galaxy with average inclination as a func- 
tion of HI velocity width. The correction is im- 
plemented only when comparing the observations 
with our model galaxies. 

2.2. Early- type galaxies 

We also include bulge-dominated early-type 
galaxies (ellipticals and lenticulars) in the LV re- 
lation, again measuring the circular velocity at 
our fiducial 10 kpc radius. The circular velocity 
in this case is used not as a measure of rotation 
but merely as a probe of the mass profile, fur- 
ther justifying the use of the term "LV relation" 
instead of TF relation. Using early-type galaxies 
allows us to probe closer to the mass regime where 
the abundance of DM halos drops exponentially 
(i.e., the knee of the velocity function), which is 
more sensitive to the cosmological model. It also 
allows for study of halo-galaxy relations without 
regard to the details of the evolution of the stellar 



populations within them. 

Because of the challenges of both observing and 
modeling early-type galaxies, so far there exists 
no comprehensive set of mass measurements for 
them akin to the spiral galaxy samples. Instead, 
we compile a set of high-quality LV estimates for 
individual galaxies from the literature. 

To provide the necessary LV data for nearby 
elliptical and lenticular galaxies, we searched the 
literature for high-quality mass measurements at 
~ 10 kpc radii. A variety of different mass trac- 
ers were used including hot X-ray gas, a cold gas 
ring, and kinematics of stars, globular clusters, 
and planetary nebulae. We required the mass 
models to incorporate spatially-resolved temper- 
ature profiles in the case of X-ray studies, and to 
take some account of orbital anisotropy effects in 
the case of dynamics. We also used only those 
cases where T^io was constrained to better than 
~ 15%. We make no pretense that this is a sys- 
tematic, unbiased, or especially accurate sample of 
early-type masses, noting simply that it is prefer- 
able to ignoring completely this class of galaxies 
which dominates the bright end of the luminosity 
function. 

As an aside, we find in comparing to cen- 
tral v elocity dispersions ag taken from Hyper- 
Leda (|Paturel et al.l 120031 ). that the scaling T4o — 
v2(To works very well on average, suggesting near- 
isothermal density profiles over a wide range of 
galaxy masses. It is far easier to measure ctq ob- 
servationally than Vio, motivating the use of the 
former as a proxy for the true Vdrc which is more 
robustly predicted by theory. The ~ 15% scat- 
ter that we find in the ao-Vio relation is relatively 
small, but it is beyond the scope of this paper to 
consider the potential systematics of using (Tq as a 
proxy. We will use Vio for the LV analysis in this 
paper. 

For the luminosities, we make use of the to- 
tal B-h&nd apparent magn itudes from the RC3 
(|de Vaucouleurs et al.lll99ll) . corrected for Galac- 
tic extinction. To correct t o '^•^r magnitudes, 
we use the filter conversions in IBlanton fc Roweia 
(|2007l ) together with (B - V) colors obtained 
from HyperLeda ( Paturel et al.l 12003). For the 
distances (required both for absolute magnitudes 
and for choosing the circular velocity measurement 
radii in kpc), we use as a first choice the estimates 



from surface brightness fiuctuations (jJensen et al 
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Fig. 1. — Properties of the early- type galaxy sample. Left panel: _B-band luminosity versus circular velocity 
measured at 10 kpc for individual galaxies. Symbols indicate the mass probe used: stellar kinematics 
[light blue stars), X-ray gas (purple crosses), planetary nebula kinematics (green squares), globular cluster 
kinematics (orange filled circle), and a cold gas ring (red open circle). Dashed lines show _B-band dynamical 
mass-to-light ratios of M/Lb = 3, 6, 12, and 24 (top to bottom) assuming all the light is contained within 10 
kpc; for comparison, early-type galaxies are expected to have stellar M/Lb ~ 2.0-2.3. Right panel: Stellar 
mass as a function of circular velocity at 10 kpc for ellipticals and SOs along with some late types shown for 
reference. 



20031 ). and otherwise the recession velocity. 

The local data for 55 individual early-type 
galaxies are presented in Figure [1] (left). Dashed 
lines show i?-band dynamical mass-to-light ra- 
tios of M/Lb = 3, 6, 12, and 24; for com- 
parison, early-type galaxies with typical colors 
(B — V ^ 0.85-0.95) are expected to have stel- 
lar M/Lb ^ 2.0-2.3 for a Chabrie r IMF (e.g., 
Fig. 18 of iBlanton fc Roweid (|2007t )V A table 
including the sources of the data is provided in 
Appendix |B] There is no obvious systematic dif- 
ference between the results from different mass 
tracers. The galaxies appear to trace a fairly tight 
LV sequence, except around the L* luminosity 
(assuming M^ ss —20.6), where there are a few 
galaxies whose circular velocities appear to be 
fairly high or low. The low-Vio galaxies include 
NGC 821 and NGC 4494, which were previously 
suggested as having a "dearth of dark matter" 



(jRomanowskv et al.ll2003l ) , and as implying a dark 
matt er "gap" with respect t o X-ray bright ellipti- 
cals ( Napolitano et al.l 120091 ). The present compi- 
lation suggests that the galaxy population in the 



local universe may fill in this gap, although further 
work will be needed to understand the scatter. 

The right panel of Figure [T] shows the relation 
between stellar mass and circular velocity at 10 
kpc for the galaxies in the early- type sample along 
with a few spirals for comparison. Stellar masses 
were obtained using equation ([8]) as explained in 
Section [^751 The ellipticals are virtually indistin- 
guishable from the SOs in the regime where they 
overlap while the spirals seem to contain slightly 
more stellar mass at the same Kirc We will dis- 
cuss this issue in more detail in Section 16.31 . 

2.3. Observational LV relation 

Figure [2] shows the combined LV relation for 
galaxies with very different morphologies: from 
Magellanic dwarfs with Vcirc ~ 50 km s^^ to gi- 
ant ellipticals with Vcirc ~ 500 km s^^. The LV 
relation is not a simple power-law. Dwarf galax- 
ies show a tendency to have lower luminosities as 
compared with a simple power-law extrapolation 
from brighter magnitudes. There is a clear sign 



of bimodality at the bright end of the LV relation 
with early type galaxies having ~ 20 — 40% larger 
circular velocities as compared with spiral galaxies 
with the same r-band luminosity (or, conversely, 
~ 1 magnitude fainter at fixed ycirc)- 
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Fig. 2. — The r-band luminosity versus circular 
velocity for several observational samples across 
different morphological types. All absolute mag- 
nitudes were transformed to the SDSS r-band at 
redshift z — 0.1. Points with error bars show the 
median and 1-a scatter. 



The Bolshoi simulation 



the 



The Bolshoi simulation was run using 
following cosmological parameters: fimattor = 
0.27, h = 0.70, (T8 = 0.82, f^bar = 0.0469, 
n = 0.95. These parameters are compatible 
with the WMAP seven-year data (WMAP7) 



(jjarosik et al.l l2010t ) and with WMAP 5 com- 
bined with Baryon Acoustic Oscillations and 



Type la Supernova data ([Hinshaw et al.l l2009t 
Komatsu etall 120091: iDunklev et al.l l2009l) . The 
Bolshoi parameters are in excellent agreement 
with the SDSS.niaxBCG+WMAP5 cosmological 
parameters ( Rozo et al.ll2009l ) and with cosmolog- 
ical parameters from WMAP 5 plus rece nt X-ray 
cluster survey results ( Klvpin et al.ll2010l ). 

It is important to appreciate that Bolshoi dif- 
fers essentially from another large, DM-only cos- 



molo gical simulation: t he Millennium simula- 
tion (|Springel et al.ll2005l MS-I). The Millennium 
simulation has been the basis for many studies 
of the distribution and statistical properties of 
dark matter halos and for semi-analytic models 
of the evolving galaxy population. This simula- 
tion and the more recent Mi llennium-II simulation 
(JBoylan-Kolchin et al.ir2009l MS-II) used the first- 
year (WMAPl) cosmological parameters, which 
are rather different from the most recent esti- 
mates. The main difference is that the Millennium 
simulations used a substantially larger amplitude 
of perturbations than Bolshoi. Formally, the value 
of (78 used in the Millennium simulations is more 
than 3cr away from the WMAP5+BA0+SN value 
and nearly 4cr away from the WMAP7-|-BAO-(-Ho 
value. However, the difference is even larger on 
galaxy scales because the Millennium simulations 
also used a larger tilt n = 1 of the power spectrum. 
The Bolshoi simulation uses a computational 
box 250 /i~i Mpc across and 2048^ « 8.6 biUion 
particles, which gives a mass resolution (one par- 
ticle mass) of TOi = 1.35 x 10^ /I'^M©. The force 
resolution (smallest cell size) is physical (proper) 
1 h~^ kpc. For comparison, the Millennium-I 
simulation had a force resolution (Plummer soft- 
ening length) 5 h~^ kpc and the Millennium- 
II simulation had 1 h^^ kpc. The Bolshoi sim- 
ulation was run with the Adaptive-Refinement- 
Tree (ART) code, which is an Adaptive-Mesh- 
Refinement (AMR) type code. A detailed de- 
script io n of the code is given in iKravtsov et al. 
( 19971): iKrav tsov (1999). We refer the reader to 



Klvpin et al.l ( 20101 ) for more details specific to the 



use of the code for the simulation. 

We use a parallel (MPI-|-OpenMP) version 
of the Bound-Density-Maxi ma (BDM) algorithm 
to id entify halos in Bolshoi (JKlvpin fc Holtzman 
19971 ). BDM does not distinguish between halos 



and subhaloqll - they are treated in the same way. 
The code locates maxima of density in the dis- 
tribution of particles, removes unbound particles, 
and provides several statistics for halos including 
virial mass and radius, as well as density profilef|j. 
We use the virial mass definition Afvir that follows 
from the top-hat model in the expanding Universe 



^A subhalo is a halo which resides within the virial radius 

of a larger halo. 
^The Bolshoi halo catalogs are publicly available at 

|http : //www . mult Idark ■ org | 



with a cosmological constant. We define the virial 
radius i?vir of halos as the radius within which the 
mean density is the virial overdensity times the 
mean universal matter density pm = ^mPcrit at 
that redshift. Thus, the virial mass is given by 



-^^vir = -^^virPm-Rvii 



(1) 



For our set of cosmological parameters, at z = 
the virial radius i?vir is defined as the radius of 
a sphere enclosing average overdensity equal to 
Avir = 360 times the average matter density. 
The overdensity limit changes with redshift and 
asymptotically goes to 178 for high z. Different 
definitions are also found in the literature. For 
example, the often used overdensity 200 relative 
to the critical density gives mass M20O) which 
for Milky- Way-mass halos is about 1.2-1.3 times 
smaller than Afvir- The exact relation depends on 
halo concentration. 

At each timestep there are about 10 million ha- 
los in Bolshoi (8.8 xlO^ at z = 0, 12.3 xlO^ at z = 
2, 4.8 xlO^ at z = 5). The halo catalogs are com- 
plete for halos with Vc\xc > 50 km s^^ (Mvir ~ 
1.5 X 10^° /i^^Mq). In order to track evolution 
of halos over time, we find and store the 50 most 
bound particles. Together with other parameters 
of the halo (coordinates, velocities, virial mass, 
and circular velocity) the information on most 
bound particles is used to identify the same ha- 
los at different moments of time. The procedure 
of halo tracking starts at z = and goes back in 
time. The final result is the history (track) of the 
major progenitor of a given halo. 

4. DM halos: definitions and characteris- 
tics 

We distinguish between two types of halos. A 
halo can be either distinct (not inside the virial ra- 
dius of a larger halo), or a subhalo if it is inside of 
a larger halo. For both distinct halos and subha- 
los, the BDM halo finder provides the maximum 
circular velocity 



Vci 



GM{< r) 



(2) 



Throughout this paper we will use term circular 
velocity to mean maximum circular velocity. 



As the main characteristic of the DM halos we 
use their circular velocity Vcirc- There are advan- 
tages to using Vcirc as compared with the virial 
mass Mvir ■ The virial mass is a well defined quan- 
tity for distinct halos, but it is ambiguous for sub- 
halos. It strongly depends on how a particular 
halo-finder code defines the truncation radius and 
removes unbound particles. It also depends on the 
distance to the center of the host halo because of 
the effects of tidal stripping. Instead, the circu- 
lar velocity is less prone to those complications. 
The main motivation for using V^irc in this work 
is that it is more closely related to the proper- 
ties of the central regions of halos and, thus, to 
galaxies hosted by those halos. For example, for 
a Milky- Way type halo the radius of the maxi- 
mum circular velocity is about 40 kpc (and Vdrc 
is nearly the same at 20 kpc), while the virial ra- 
dius is about 200 kpc. In addition, the virial mass 
of a DM halo is not an easily observable quantity 
and this further limits its use for comparison of 
simulations with observations. 

Tidal stripping can lead to significant mass loss 
in the periphery of subhalos. The net effect at 
redshift zero of the complex interactions that each 
halo undergoes is a decrease in the maximum cir- 
cular velocity compared to its peak value over 
the entire history of the halo. The galaxy re- 
siding in the central region of the halo should 
not experience much stripping and should preserve 
most of its mass inside the optical radius (e.g., 
Nagai fc Kravtsovll2005HConrov et allbOOGl ). Fol- 



lowing this argument, the initial total mass distri- 
bution and rotation profile of the halo are frozen at 
the moment before the halo is accreted and starts 
to experience stripping. We refer to this circu- 
lar velocity as Vacc- In practice we find the peak 
circular velocity of the halo over its entire history. 
Further det ails on the halo track ing procedure can 
be found in lKlvpin eraP (|2010l ). 



5. Connecting galaxies and DM halos 

To investigate the statistics of galaxies and their 
relation to host DM halos as predicted by the 
ACDM model using HAM, we obtained the prop- 
erties of our model galaxies using the following 
procedure: 

1. Using the merger tree of each DM halo and 
subhalo, obtain Vacc = the peak value of 



the circular velocity over the history of the 
halo (this is typically the circular velocity of 
the halo when it is first accreted). Perform 
abundance matching of the velocity function 
of the halos to the LF of galaxies to obtain 
the luminosity of each model galaxy. 

2. Perform abundance matching of the veloc- 
ity function to the stellar mass function of 
galaxies to obtain the stellar mass of each 
model galaxy. 

3. Use the observed gas-to-stellar mass ratio as 
a function of stellar mass to assign cold gas 
masses to our model galaxies. The stellar 
mass added to the cold gas mass becomes 
the cold baryonic mass. 

4. Using the density profiles of the DM halos, 
obtain the circular velocity at 10 kpc (Vio) 
from the center of each halo. To do this, cal- 
culate the dark-matter-only contribution by 
multiplying the DM mass profile, obtained 
directly from the simulation, by the factor 
(1 ^ /bar); where /bar IS the cosmological 
fraction of baryonqj- Then take the total 
cold baryon contribution from step 3 and as- 
sume it is enclosed within a radius of 10 kpc. 
Adding the two contributions gives the total 
mass required to calculate Vlo. 

5. Implement the correction to Vio due to the 
adiabatic contraction of the DM halos due 
to the infall of the cold baryon component 
to the center. 

We now explain each of the above 5 steps in 
detail. 

Using the key assumption that halos with 
deeper potential wells become sites where more 
baryons can gather to form larger and more lumi- 
nous galaxies, we ranked our halos and subhalos 
using their Vcirc , and starting from the bright end 
of the r-band LF, assigned luminosities to each 
according to thei r space density usin g the pre- 
scription found in IConrov et al.l ( 20061 ) . In other 
words, we found the unique one-to-one correspon- 
dence that would match the halo velocity function 



with the luminosity function of observed galaxies. 
Of course, this is a simplifying approximation. It 
does not discriminate between blue (star-forming) 
and red (quenched) galaxies, for example. 

In this paper we use the Sch echter fit to the r- 

band g alaxy LF measurement of M ontero-Dorta fc Prada 
(I2OO9I) obtained from the SDSS Data Release 6 



(DR6) galaxy samplqj. The fit is characterized by 
the parameters: $* = 0.0078, M*-^^ - 5 log ft. = 
—20.83, and a = —1.24. Since the median redshift 
of the SDSS DR6 sample is « 0.1 (JBlanton et al. 



2003af ). all our subsequent results will be shown 



in r-band magnitudes. 

As an alternative, we also consider a LF with a 
steepe r slope at low-luminosities. 



Blanton et al 



(|2005l ) obtained the SDSS LF including dwarfs 
as faint as Mr = — 12 and investigated surface 
brightness completeness at the faint end of the 
distribution. Their steeper value of the faint-end 
slope was obtained by weighting the abundance 
of each galaxy by its estimated surface brightness 
completeness. To quantify the effect of including 
low surface brightness galaxies in our model (those 
with Petrosian half-light r-band surface brightness 
greater than 24.0 mag/arcsec^), we increased the 
abundance of bright gala xies in the Blanton et al. 
( 2005 ) LF to match the iMontero-Dorta fc Prada 
(J2OO9I ) LF at the bright end while keeping the 
steep faint-en d slop e (a = —1.34) measured by 



Blanton et al 



(|2005l ). The modified LF produces 
60% more galaxies at Mr ~ — 16 and over a factor 
of 2.5 more galaxies at Mr ~ —13. Using this LF 
to perform the abundance matching increases the 
luminosity of galaxies assigned to small DM halos, 
steepening the faint end of the LV relation. 

It is important to note that we assume that 
each dark matter halo or subhalo must contain 
a galaxy with a detectable luminous component 
(for the SDSS r-band this requires galaxies to be 
detectable in visible wavelengths) and this com- 
ponent must evolve in a way that guarantees its 
detectability at z = 0. Since the effective vol- 
ume surveyed by SDSS DR6 at z < 0.3 is com- 
parable to the volume of the Bolshoi simulation, 
we expect the statistics of the halo population to 
be comparable to those of the observed galaxies 
all the way up to the large mass/luminosity tail 



•^ Recall that the Bolshoi simulation was run for a dissi- 
pationless cosmic density Cm = f^dm "I" f^^ar = 0.27 = 

f^dm(l + /bar)- 



*To avoid aperture corrections when comparing to other 
data we use model magnitudes instead of Petrosian values. 



of the distributions. Even though Bolshoi con- 
tains a factor of ^ 1.8 more objects th an the sam- 
ple of iMontero-Dorta fc Pradal ( 20091 ) . abundance 
matching is mostly insensitive to uncertainties in 
the high-luminosity tail of the LF. 

5.1. The circular velocity of galaxies inside 
halos including cold baryons 

The next step is to separate the DM and baryon 
components in each halo and allow the baryons to 
dissipatively sink to the centers of the DM halos. 
We assume for simplicity that there is a radius at 
which we could consider most of the cold baryons 
to be enclosed, with only dark matter present be- 
yond that point. 

This cold baryon component has been observed 
to comprise only a small fraction of the cos- 
mic abundance of baryons; in other words, the 
cold baryon fraction /bar = Mbar/(AfDM + Mha.r 
in galaxies is much lower (JFukugita et al.l Il99 



ization of the observed atomic gas mass fraction 



as a function of stellar mass from iBaldrv et al 



[Fukugita & Peebles 



0.17 (Komat su et al 



200 



200' 



than Sibar/' 'matter = 

We resort to the ob- 
servations and use the galaxy stellar mass func- 
tion o btained from the SDSS DR7 bv lLi fc Whit3 
(12009), who employ e stima tes of stellar masses 
bv iBlanton fc RoweisI (|2007l ) obtained using five- 
band S DSS pho t ometr y assuming the universal 
IMF of lChabrieii ( 20031 ). These masses are consis- 
tent with those estimate d using single-colo r and 
spectroscopic techniques (JLi fc Whitell2009() . 

Using the same procedure described above for 
the luminosity function, we abundance-matched 
the halos in Bolshoi to the galaxies in the 
SDSS DR7 starting from the high stellar mass 
end until reaching our completeness limit at 
T4irc — 50 km s~^, obtaining stellar masses for 
each galaxy. Strictly speaking, this procedure re- 
sults in a one-to-one relation between circular ve- 
locity and stellar mass-to-light ratio which should 
only be interpreted as the average of a population 
of galaxies with a given Vdrc- The scatter (or bi- 
modality) in the mass-to-light ratio as a function 
of circular velocity could be measured from obser- 
vations and included in the assignment but would 
not change our results significantly. 

Since dwarfs can have most of their cold 
baryons in the ga s phase instead of in stars 
(jBaldrv et al.ll2008l ). we calculated for each model 
galaxy the total cold gas mass using a parameter- 



(J2008I ) (their equation (9), shown as a dashed line 
in their Figure 11). This includes the total cold 
atomic gas found in the disk only. The gas-to- 
stellar mass ratio /gas depends on stellar mass 
and it is the largest for dwarfs. For example, 
/gas ~ 4 — 5 for galaxies with M* = 10* M©. 
It declines to /gas ~ 0.25(0.1) for galaxies with 
M^ = 10i°(10") Mq. It should be even smaller 
for ellipticals and SOs. It should be noted that, 
when it comes to dynamical corrections to Vcirc, 
the gas plays a minor role. It only becomes impor- 
tant when we consider the baryonic Tully-Fisher 
relation. 

Lastly, we add the stellar and cold gas masses 
for each model galaxy and obtain the correction 
to the circular velocity of the pure DM halo at a 
radius enclosing the cold baryonic mass. We set 
this value to 10 kpc for all the halos in our sample. 
In the case of dwarf galaxies this should be a good 
approximation to the maximum circular velocity 
since their rotation curves rise much more slowly 
and in some cases they peak beyond the optical 
radius (jCourteaul 119971) . Our assumption allows 
us to include the peak of the rotation curve for 
most of these galaxies. In the case that the peak 
is located well within 10 kpc the correction would 
be almost negligible since we would be still mea- 
suring rotation in the flat regime. Additionally, 
truncating the cold baryons at a radius of 10 kpc 
allows us to directly calculate the correction to the 
circular velocity at that radius without having to 
resort to more complicated assumptions about the 
distribution of baryonic matter in galaxies, i.e., ex- 
ponential lengths and Sersic indices of disks and 
bulges as well as extended gas and stellar halos. 

To obtain the circular velocity measured at 
10 kpc (Vio) for the Bolshoi DM halos, we need 
to use dark matter profiles and find the dark mat- 
ter mass inside a 10 kpc radius. To do this, we 
could use the individual profile of each halo. How- 
ever, once we select halos with a given maximum 
circular velocity, individual halo-to-halo variations 
are small at 10 kpc (the situation is different if 
we select halos using virial mass, which results in 
large deviations in concentration producing large 
variations in Vio). This is why instead of individ- 
ual profiles we construct average profiles for halos 
within a narrow range A logj^Q V ~ 0.05 of maxi- 



10 



mum circular velocity. 

We first bin and average the circular velocity 
profiles of the distinct halos found by the BDM 
code. These profiles are calculated for each halo 
(including unbound particles) in logarithmic ra- 
dial bins in units of i?vir- Using distinct halos is 
convenient because it gives us density profiles that 
are less affected by interactions than those of sub- 
halos. For the inner profiles of subhalos the effect 
is relatively small because the stripping happens 
preferentially at the outer radii. Using the av- 
eraged binned circular velocity profiles we obtain 
the velocity at 10 kpc. Within about 1.2% of the 
virial radius, discreteness effects render the pro- 
files unreliable and we use instead extrapolation 
with the shape of a simple power-law in radius. 
For halos with Vchc < 100 km s~^ the full ex- 
tent of the profiles suffer from measurement noise 
which we avoid by extrapolating from the profile 
of halos with ~^ 100 km s~^. Figure [3] shows the 
obtained median relation between the maximum 
circular velocity V^nax and V^io for the Bolshoi DM 
halos without inclusion of the cold baryons. 

The estimates of the relations obtained when a 
parametric form of the density pro file is used are 



also shown for the case of the NFW ( Navarro et al 



19971) 



p(r) 



4ps 



(3) 



and the Einasto ("E inastd Il965t iGraham et al.l 
20061 ) universal profiles 



p{r) = ps exp {-2n [(r/r^)!/" - l] } ; (4) 

where rs is the radius at which the logarithmic 
slope of the density profile is —2. Following 



Graham et al.l (|2006f ). we use n — 6.0. The con- 
centration parameter defined for both models as 
c = Ryir/rs is given by th e relations obtaine d in 



Pap er I for distinct ha los (JKlvpin et al.ll20ld see 
also IPrada et al.1l201ll) 



9.60 



1012 h^'^Mr. 



-0.075 



(5) 



and 

M ■ — 



Vr; 



2.8 X 10-2 km s" 



3.16 



h-'Mr^. (6) 



Note that in Figure [3] we use total dynamical 
masses and do not account for the condensation 
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Fig. 3. — Relation between the maximum circular 
velocity (Vmax) and the circular velocity measured 
at 10 kpc (Vio) for the dark matter halo only (ex- 
cluding the cold baryonic component). The solid 
curve shows the binned median values of the Bol- 
shoi DM halo sample. The other curves show the 
relation obtained assuming the NFW (dashed) and 
the Einasto (dot-dashed) profiles with the halo 
concentration given by eq. ([S|). 

of baryons. For Vdrc = 100 — 450 km s^^ the 
rotation (or density) profiles of the Bolshoi sim- 
ulation are extremely well approximated by the 
Einasto parameterization, whereas NFW underes- 
timates T4o by almost 20% at 450 km s~^. Fol- 



lowi ng the conclu s ions of iNavarro et al.l (J2004l ) 



and lGraham et al.l (|2006[ ). we assume the Einasto 
profile when extrapolating the inner parts of the 
largest {Vdrc > 450 km s^^) halos. 

5.2. Adiabatic contraction of DM halos 

Dissipation allows the baryons to condense into 
galaxies in the central regions of DM halos drag- 
ging the surrounding dark matter into a new 
more concentrated equilibrium configuration. If 
the density of the DM halo increases consider- 
ably within the extent of the disk, the peak cir- 
cular velocity could be much larger than our pre- 
vious estimates. There are different approxima- 
tions for the adiabatic compression of the dark 
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Fig. 4. — Comparison of the observed Luminosity- Velocity relation with the predictions of the ACDM model 
using halo abundance matching. The solid curve shows the median values of °-^r-band luminosity vs. circular 
velocity for the model galaxy sample. The shaded regions enclose 68% and 95% of the model galaxies in bins 
of luminosity. The circular velocity for each model galaxy is based on the peak circular velocity of its host 
halo over its entire history, measured at a distance of 10 kpc from the center including the cold baryonic mass 
and the standard correction due to adiabatic halo contraction. The dashed curve shows results for a steeper 
(a = —1.34) slope of the LF. The dot-dashed curve shows predictions after adding the cold baryon mass 
but without adiabatic halo contraction. Points show median and 1-cr scatter of representative observational 
samples. 



matter. iBlumenthal et al.l ( 19861 ) provide a simple 
analytical expression, which is known to overpre- 
dict the effec t . Th e approximation proposed by 
ICnedin et al.l ()2004l ) predicts significantly smaller 
increase in the density of the dark matter. More 
recen t simulations indicate even smaher co mpres- 
sion (JTissera et al.ll2009HDu"ffv et al.ll2010l) . How- 
ever, at a 10 kpc radius the dark matter con- 
tributes a relatively large fraction of mass even for 
large galaxies. As a result, the difference between 
the str ong compression model of lBlumenthal et al 



use the standard adiabatic contraction 



( 19861 ) and no-compression is only 10 — 20% in ve- 



locity. 

We 

(AC) model of lBlumenthal et al.l (jl986[ ) to bracket 
the possible effect. We thus assume that following 
the condensation of the baryons, the dark matter 
particles adjust the radius of their orbits while 
conserving angular momentum in the process. We 
solve the equation 

Mtot(n)n = [MDM(n)(l - /bar) + Mbar(rf)]rf, 

(7) 
where ri — 10 kpc, Afbai(n) is the total bary- 
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onic mass assigned to each halo and /bar = 
f^bar/^mattor IS the Universal fraction of baryons. 
We solve equation ^ for n and then add the 
dark matter mass MDM(n)(l — /bar) to the mass 
of cold baryons to find the circular velocity. Note 
that this implies that only cold baryons (i.e., stars 
and cold gas) are left in the central regions of 
the galaxy, while the remaining hot baryons are 
at larger radii. As expected, only the halos that 
are dominated by baryons at their centers suffer 
a significant increase in their measured circular 
velocities due to the increase in concentration of 
dark matter as result of adiabatic contraction. 

6. Results 

6.1. The Luminosity- Velocity relation 

Figure |4] shows the predicted LV relation for 
galaxies in the ACDM model obtained using halo 
abundance matching. We binned together the ma- 
jor observational samples described in Section [5] 
and include them for comparison. The inter- 
nal extinction de-correction for late-type galax- 
ies described in Section 12.11 was implemented for 
a fair comparison with the models. The full 
curve in this plot shows r esults obtained using the 
Montero-Dorta fc Pradal ([2009) LF of galaxies in 
the SDSS DR6 sar nple and uses the adiabat ic con- 



(|l98fih . The 



traction model of iBlumenthal et al, 
1- and 2-(j width of the distribution of model 
galaxies in bins of luminosity is represented by 
the shaded regions (details about the model used 
to add scatter can be found in Section 16.2. 4p . 
Predictions without adiabatic contraction (with 
the cold baryon contribution added in quadra- 
ture) are shown as the dot-dashed curve. The 
dashed curve shows the effect of a steeper slope 
at the faint end of the LF that accounts for poten- 
tial surface brightness incompleteness. (For de- 
tails see Section [5|). Although dwarf galaxies with 
T4irc < 80 km s ~^ seem to agree better with a 
model using the iMontero-Dorta fc Pradal (|2009l ) 
luminosity function, the scatter of the observed 
dwarf LV relation is so large that both LFs used 
in conjunction with the abundance of DM halos 
produce results that are consistent with the avail- 
able data. 

One may note that the AC model misses late- 
type points with ycirc = 150 — 250 km s^^ and the 
no- AC model practically fits most of the late-type 



measurements. This should not be interpreted as 
either an advantage for the no- AC model or a dis- 
advantage for the AC model. Our predictions ap- 
ply to the average population across all types of 
galaxies. Because of the dichotomy of the LV dia- 
gram, a model that goes through either early-type 
galaxies or through late-type galaxies is an incor- 
rect model. The correct answer should be a model 
which tends to be close to spirals at low lumi- 
nosities (where spirals dominate the statistics) and 
gradually slides towards the ellipticals at the high- 
luminosity tail where they dominate. It seems that 
the AC model does exactly that, but it may over- 
predict the circular velocities at the very bright 
tail of the LV relation^ 

As demonstrated in Section fG. 2. 41 our stochas- 
tic HAM model accounts for galaxies that reside 
in halos with both smaller and larger Vcirc than 
the average. For example, since the most massive 
spiral galaxies comprise a very small percentage 
of the galaxy population with V^irc > 250 km s~^ 
(less than 10%), in our model they are assigned 
to the low-V^irc wing of the distribution shown in 
Figure |4l Hence, even though the model makes 
predictions for the average galaxy population, it 
also accounts for the morphological bimodality ob- 
served in the LV relation. 

One also should not overestimate the quality 
of the observations. The fact that in Figure |4] the 
brightest spirals with Af^ — 5 log^^Q < —21 are more 
than 2-a away from the mean of the models is 
not a problem because of the size of the uncer- 
tainties in the observational data. For instance, 
there is a systematic ~ 10 % velocity off s et be- 
twe en the measurements o f iPizagno et al.l (|200lT ) 
and ISpringob et al.l (|2007l ). which seems to point 
to the current uncertainties in the LV relation. 

Considering the current level of the uncertain- 
ties, our model galaxies show remarkable agree- 
ment with observations spanning an order of mag- 
nitude in circular velocity (or, equivalently, three 
orders of magnitude in halo mass) and more than 
three orders of magnitude in luminosity. For 
galaxies above 200 km s~^, our model produces 
results that agree extremely well with the ob- 
served luminosities of early-types (Es and SOs). 



Schulz et ai] | |2010|) . iNapolitano et al.1 | |2010D and 
Tortora et all l|201Cl ) find observational evidence for 
AC in elliptical galaxies. 
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Given how simple the prescription is, it is per- 
haps surprising how closely we can reproduce the 
properties of observed galaxies. For galaxies with 
T4irc = 100 - 200 km s~\ DM halos without any 
corrections already occupy the region expected 
for galaxies. The dynamical corrections improve 
the fits, but it is important to emphasize that 
the abundance matching method yields the cor- 
rect normalization of the LV relation regardless of 
the details of the corrections we implement. An- 
other feature of the relation, its steepening below 
100 km s~^, is caused by the the shallow faint- 
end slope of the luminosity function. Although 
our model galaxies in this regime are slightly un- 
dcrluminous as compared with a simple power- 
law extrapolation from brighter galaxies, observed 
dwarfs seem to predict a deviation from a power- 
law TF in the same general direction. 

In the way it was constructed, our model galaxy 
sample does not include uncertainties in either the 
halo velocity function or in the galaxy luminosity 
function. This produces an LV relation with no 
scatter. Section 16.2.41 examines the effects of in- 
cluding scatter in the halo matching procedure. 

We now discuss in greater detail the results of 
the individual steps explained in Section [5] 

6.2. The Luminosity- Velocity relation: de- 
tailed analysis of model components 

6.2.1. Measuring circular velocity at 10 kpc 

Observations do not always provide measure- 
ments of the circular velocity at 10 kpc. This is 
especially true for dwarf galaxies where the last 
measured point of the rotation curve can be at 3- 
5 kpc. What are the errors associated with using 
measurements at different radii? Figure [5] presents 
three typical examples of the circular velocities ex- 
pected for galaxies with vastly different masses. 

The top panel shows a model of a giant ellip- 
tical galaxy with 1.5 x 10^^ Mq of stellar mass 
distributed according to a R^^^ law with a half- 
light radius of 5.5 kpc. The stellar component is 
embedded in a dark matter halo with virial mass 
Mvir = 10"'^^ Mq and median concentration c = 7. 
The maximum circular velocity (310 km s^^) of 
the halo is reached at 160 kpc. The middle 
panel shows a Milky Way-size model with max- 
imum circular velocity 190 km s~^, virial mass 



c = 9 for its mass. The cold baryonic component 
consists of a Hernquist bulge (Mbuigc = 10^" Mq, 
half-mass radius -Rbuigo = 1 kpc) and an expo- 
nential disk (Mdisk = 5 X 10^° Mq, scale radius 
-Rdisk = 3.5 kpc). The bottom panel shows a dwarf 
galaxy model with Mvir = 7 x 10^*^ Mq, c = 12, 
^lax — 65 km s~^. Its cold baryons have two ex- 
ponential disks: one for stars and and one for cold 
gas with a mass ratio of 1:4 and radii i?, = 1.5 kpc 
and i?gas — 3.0 kpc. The total mass in cold 
baryons is Afbar = 3 x 10^ Mq. When we include 
baryons, we assume that most of them were blown 
out from the models and the only baryons left are 
in the form of stars and cold gas. As in Bolshoi, 
the "DM" circular velocities in Figure [S] include 
a cosmological amount of baryons that traces the 
distribution of the dark matter. This contribution 
is removed from the mass profiles before adding 
the cold baryons. In all three cases we use the 
Einasto dark matter profiles (equation (U)) with 
n = 6. For the models labeled "DM+Baryons" 
at each radius we simply add the mass of cold 
baryons and the mass of dark matter. The mod- 
els termed "DM-|-Baryons-|-AC" include adiabatic 
compressio n of dark matter a- c cordin g to the pre- 
scription of iBlumenthal et al.l (|l986l ). 

After adding the cold baryons the circular ve- 
locity profiles become rather flat in the inner 
5—10 kpc regions of the galaxies implying that 
measurements of circular velocities anywhere in 
this region are accurate enough to provide the 
value of the circular velocity at 10 kpc. 

There are some caveats in choosing 10 kpc 
as a fiducial radius for either extremely massive 
spirals or giant ellipticals. Considering the cor- 
relation between central surface brightness an d 
disk scale-length found by ICourteau et al.l (|2007| ) , 
the most luminous disks appear to have scale- 
lengths as large as 15 kpc, whereas according to 
Courteaul (|l997l ) their rotation curves peak at 



A/v 



1.7 X 10 Mq, and median concentration 



about 1 scale-length. Hence, we may be underes- 
timating the maximum observed rotation velocity 
of these galaxies in our sample. We may also over- 
estimate circular velocities when we assume that 
most of the cold baryon mass is inside 10 kpc ra- 
dius. In principle, some corrections can be applied 
to compensate this effect. However, our estimates 
show that at most this is a ^ 20% effect for spirals 
and somewhat smaller for ellipticals because they 
are more compact for the same luminosity. Con- 
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Fig. 5. — Effect of cold baryons on circular ve- 
locity profiles for three characteristic models of 
galaxies with virial masses 10^^ M© (top), 1.7 x 
10^^ M0 (middle), and 7 x 10^° Mq (bottom). 
The "DM" curves include a cosmological fraction 
of baryons that trace the dark matter distribu- 
tion. The cold baryon mass is added to the true 
dark matter mass in calculating the circular ve- 
locity ("DM+Baryons"). The effect of adiabatic 
compression of the dark matter is included in the 
models named "DM-|-Baryons-|-AC" . After adding 
the cold baryons the circular velocities are rather 
flat in the inner 5 — 10 kpc regions. 



sidering existing uncertainties and complexities of 
implementing the correction, we decided not to use 
them. 

6.2.2. Dark matter profiles 

To illustrate the effect of tidal stripping. Fig- 
ure [6] shows the results (dashed curve) obtained 
when the luminosity assignment is performed us- 
ing the peak historical value of each halo's cir- 
cular velocity (Vacc) as compared with the circu- 
lar velocity at 2; = (dotted curve). The reason 
why the dashed curve is rightwards of the dot- 
ted one is that for subhalos the circular velocity 
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Fig. 6. — Comparison between different effects on 
the measured circular velocities of model galax- 
ies without corrections for cold baryons or adi- 
abatic contraction. The dotted curve shows the 
median ^'^r-band luminosity vs. circular velocity 
of the model galaxies that results from abundance 
matching using the maximum circular velocity of 
each DM halo at z = 0. The dashed line shows the 
effect of using the peak value of the maximum cir- 
cular velocity over the history of each halo (Vacc)- 
The solid curve shows the result of measuring 14cc 
for each halo at 10 kpc from the center. This af- 
fects intermediate and large halos the most since 
their circular velocity profile is still rising at that 
distance. All the curves include a cosmic baryon 
contribution that traces the dark matter. 



estimated at z = is smaller than its value at 
accretion Vacc- If we compare luminosities at the 
same circular velocity, the differences can be sub- 
stantial: almost one magnitude for galaxies with 
Kirc = 50—60 km s~^due to the steep slope of the 
LV relation for dwarfs. In terms of velocities, the 
differences are much smaller. Neglecting the ef- 
fects of stripping in the assignment scheme affects 
mostly dwarf galaxies, overestimating their circu- 
lar velocities by a maximum of ^ 20%. For larger 
galaxies the effect decreases to less than 5%. This 
is due to the fact that stripping only affects subha- 
los and they comprise only a minority (about 20%) 
of the total halo population. In addition, only the 
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small number of subhalos which orbit close to their 
host halo's center get significantly stripped and 
experience a substantial decline in their circular 
velocity. 

Comparison of LV relations constructed using 
T4cc one with the maximum the circular veloc- 
ity (the dashed curve in Figure |6]) and another 
with velocities T4cc estimated at 10 kpc (the full 
curve) indicates that this affects the largest halos 
the most. For example, the Vlo velocity is almost a 
factor of two smaller than Vmux for the group-size 
halos presented in the plot. Taking the circular ve- 
locity at 10 kpc also makes the LV relation much 
less curved as compared with the maximum ve- 
locity. Below ~ 80 km s~^ the maximum circular 
velocity of the DM halo happens near or within 
10 kpc, which explains why the curve does not 
shift in this regime. 

6.2.3. The effects of cold baryons and adiabatic 
compression 

Figure [7] shows how cold baryons change the 
circular velocity at a 10 kpc radius. Here we use 
two extreme approximations that bracket the ef- 
fect. The first approximation assumes that there 
is no change in the distribution of the dark mat- 
ter. All simulations so far indicate that there is 
some compression. Hence, the no compression ap- 
proximation definitely underpredicts the circular 
velocity V^iq. The second appro ximation uses the 
adiaba tic compression model of iBlumenthal et al 



(|l986l ) which produces the largest increase in the 
density of the dark matter. (The full and dashed 
curves were already shown in Figure U) . There 
are some differences between the LV relations pre- 
dicted by those approximations. However, the 
largest effect is just adding the cold baryons in 
quadrature to the circular velocity of the dark 
matter. Adiabatic compression increases the cir- 
cular velocity even further, but the effect is rela- 
tively minor because the fraction of cold baryons 
gets progressively smaller for larger galaxies. The 
amount of cold baryons used for the models is cru- 
cial for this test. As we discuss in Section [Ql the 
abundance matching predicts relatively small cold 
baryon masses for dwarfs and giants, and this is 
why the adiabatic compression in Figures [5] and 
[7] is 10-20% at the most. Again, dwarfs below 
~ 80 km s~^ are insensitive to the cold baryon 
presence. Here the full curve (DM+baryons) is 



slightly below the Vdrc of the dark matter curve 
because the latter includes a cosmological fraction 
of baryons which trace the DM, most of which are 
assumed to be blown out from galaxies (see Sec- 
tion O. 

Figure |8] shows the LV relation that we would 
obtain by assuming instead that half of all baryons 
within the virial radius or equivalently 8% of the 
virial mass are retained and are used to build the 
central galaxy (while luminosity is not affected). 
Both the shape and the normalization of the LV 
relation are incorrect, with the circular velocities 
systematically larger than observations by up to 
50%. Clearly, it is difficult to obtain the observed 
LV relation assuming a constant cold baryon frac- 
tion in the framework of the ACDM cosmology. 

6.2.4. The effects of including scatter in luminos- 
ity at fixed halo circular velocity 

So far, the abundance matching procedure we 
have used assumed a monotonic one-to-one rela- 
tion between halo circular velocity and galaxy lu- 
minosity or stellar mass. This assumption pro- 
duces average relations that can be compared 
with the medians of the observations. As shown 
by previous studies, a more detailed treatment 
of the scatter between halo and galaxy proper- 
ties may yield average relations of the brightest 
galaxies that deviate significantl y from the case 



with no scatter. For instance, Tasitsiomi et al 



(|2004l ) showed that iteratively introducing log- 
normal scatter of width 1.5 mags in the assign- 
ment of luminosities to DM halos produces an av- 
erage TF relation with massive galaxies that are 
brighter by about one magnitude compared to the 
monotoni c assignment. By treating the scatter an- 
alytically, iBehroozi et al.l ( 20101 ) found that per- 
forming halo abundance matching using their pre- 
ferred value of 0.16 dex of log-normal scatter re- 
duces the average stellar mass assigned to massive 
halos with total masses > 10^'^ Mq by up to 70% 
(when binning using virial mass) but does not af- 
fect less massive galaxies below the knee of the 
stellar mass function. 

Appendix |^ gives a detailed description of the 
method we employ to introduce scatter in our 
model. In short, we obtain luminosities for each of 
the galaxies in our sample by stochastically scat- 
tering the values obtained in the monotonic as- 
signment while forcing the preservation of the ob- 
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Fig. 7. — Effects of baryons on the LV relation. 
The dashed curve shows the circular velocity Vio 
after adding the baryon mass at the center of each 
halo without any adiabatic contraction of the dark 
matter. The solid curve shows the result of imple- 
menting the correcti on due to the adiabatic con- 
traction of the halos ( Blumenthal et al.lll986[) . For 
reference, the dotted line shows the circular ve- 
locity of the DM measured at 10 kpc (assuming 
the baryons trace the DM distribution; same as 
solid curve in Figure |6l). Baryons have little effect 
on dwarfs (V"circ < 100 km s~^) since dwarfs are 
dominated by DM beyond a few kiloparsecs. Just 
adding the baryons in quadrature has the great- 
est effect with the adiabatic compression giving a 
10-15% correction for bright galaxies. 



served luminosity function. When scattering the 
values of luminosity we do not constrain the shape 
of the probability distribution (e.g., log- normal) or 
require its width to be constant for all circular ve- 
locities. This is well justified since the shape of 
the intrinsic scatter is more difficult to constrain 
observationally (e.g., due to observational system- 
atics) . 

One parameter that our model does not cur- 
rently predict is the width of the probability dis- 
tribution of luminosity at a fixed halo circular ve- 
locity. This scatter originates from three main 
sources. The first is the observational error in the 



Fig. 8. — Effect of excessive cold baryon mass. We 
assume that half of the universal baryon fraction 
within each halo forms its galaxy. Median values 
(solid curve) of "^'^r-band luminosity vs. circular 
velocity of our model galaxies measured at 10 kpc 
from the center and including the correction due 
to adiabatic halo contraction. For comparison, the 
dotted line and the symbols reproduce the model 
and the observational data shown in Figure |H The 
model with 50% cold baryon fraction systemati- 
cally predicts galaxies that are too concentrated 
and fails to fit the observations. 



determination of the true luminosities of galaxies. 
Since we use the LF from the SDSS spectroscopic 
sample, these are the sum of the errors in pho- 
tometry plus the errors in the distances obtained 
from spectroscopy. At the mean redshift of the 
sample these combined errors are expected to be 
typically much less than 0.1 mags. The second 
source of scatter is the one present in the intrinsic 
relation between halo circular velocity and galaxy 
luminosity due to variatio n in the physica l pro- 
cesses of galaxy formation. IVerheijenI (J200ll ) stud- 
ied the nature of the intrinsic scatter in the Tully- 
Fisher relation from HI observations and obtained 
a value of aff — 0.38 in the i?-band. This value 
is consistent with the distribution in the obser- 
vational samples used in this work. Lastly, since 
we assign luminosities that are uncorrected for the 
inclination of disk galaxies, we also need to in- 
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elude the scatter that results from the distribution 
of dust extinction corrections observed in SDSS. 



Mailer et all (2009) find a fit to this distribution 
as a weak function of r-band luminosity and disk 
scale length. We adopt the value af}*" = 0.28 they 
use for a galaxy with Mk — —20. We neglect 
the errors due to photometry and distances and 
add the remaining two contributions in quadrature 
to obtain ctm,. ~ 0.5, which we use to introduce 
scatter to the model galaxies below the knee of 
the LF. Above this luminosity, where early-types 
dominate, we assume that the lack of significant 
internal extinction slightly reduces the scatter to 
~0.3. 

Figure ini shows the luminosity-binned distribu- 
tion of model galaxies in the LV relation obtained 
from the stochastic HAM scheme and compares 
it to the monotonic assignment discussed in Sec- 
tion 16.11 Since we are left with a choice regarding 
which quantity to average over, we choose to bin 
in r-band magnitude to be consistent with the bin- 
ning of the observations. The mean relation is al- 
most identical to the case with no scatter for galax- 
ies below 200 km s~^, while it becomes brighter 
by up to 0.3 magnitudes for more massive galax- 
ies. Galaxies below L* show a distribution of lumi- 
nosities that is close to gaussian as far as 2-cr away 
from the mean but has a slightly longer bright tail. 
Galaxies brighter than L* show a trend of narrow- 
ing of the distribution as well as a skewness that 
reduces the number of upscattered galaxies with 
increasing luminosity. 

To check the consistency of our approach we 
also calculated the average LV relation of our 
model galaxies obtain ed using the deconvol ution 
method described in iBehroozi et al.l ( 20101 ) and 
log-normal scatter. This procedure yields a devi- 
ation of the mean relation for Vcirc > 200 km s~^ 
towards higher luminosities that depends on the 
assumed width of the scatter. Figure [17] in Ap- 
pendix |^ shows this effect for constant CMr = 0.5 
(left panel) and dMr- decreasing from 0.5 to 0.3 
past L* (right panel). Even with a variable 
width that mimics our approach, the luminosity- 
binned spread obtai ned with the method of 
Behroozi et al.l ( 20101 ) is unrealistically large at 
the bright end. 

The differences between the results of the two 
methods actually reside in the assumptions about 
the shape of the spread. Since our stochastic as- 




Fig. 9. — The LV relation of the Bolshoi galaxies 
obtained using the stochastic assignment method 
described in Appendix El to add scatter. The solid 
(dotted) lines show the median (average) of the 
circular velocity (Vio) in bins of r-band luminos- 
ity. The shaded areas encompass 68% and 95% of 
the galaxies in each luminosity bin. The dashed 
line shows the result of monotonic assignment with 
no scatter. The small, < 0.3 mag decrease of the 
average luminosity of the brightest galaxies is op- 
posite in sign to the one obtained by binning in 
Vio- This is merely a result of binning bias due to 
the steepness of the velocity function. 



signment scheme does not constrain the scatter 
distribution to be log-normal and centered on the 
monotonic relation, the resulting skewness beyond 
L* allows it to preserve the median LV relation 
of the scatterless sample. In addition, without 
a skewed distribution it is extremely difficult to 
obtain a narrower distribution of galaxies at the 
bright end of the LV relation. 

From this analysis we conclude that the intro- 
duction of scatter in luminosity at a given halo 
circular velocity yields a median relation at the 
bright end that is sensitive to the shape and width 
of the probability distribution function used. The 
median LV relation is thus robust to uncertainties 
in the nature of the scatter below the shoulder of 
the LF, allowing for a direct comparison with ob- 
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servations. We prefer our scatter method for two 
reasons. First, it exactly preserves the luminos- 
ity function while the deconvolution method only 
does so approximately. Second, it naturally pro- 
duces an observed luminosity-binned distribution 
that becomes narrower for the brightest galaxies, 
in agreement with that expected from observa- 
tions. 

6.3. Baryon fraction and the baryonic 
Tully-Fisher relation 

For the LV relation the cold baryons played 
an ancillary role: they provided a correction to 
the circular velocity at 10 kpc. The correction is 
small for galaxies below 100 km s~^. For large 
galaxies the cold baryon contribution increases 
and typically is about half of the mass within 
10 kpc. Regardless of their role in the LV rela- 
tion, baryons are one of the prime subjects for 
the theory of galaxy formation. Unfortunately, ac- 
curate measurements of baryonic masses are also 
prone to some uncertainties. Dynamical mea- 
surements of the baryonic component are difh- 
cult because of dark rnatter- baryon degeneracies 
(e.g.. lDehnen fc Binnevlll998[ ). In other words, the 
baryon mass depends on what is assumed about 
the dark matter. Population synthesis provides 
an independent estimate of the stellar mass, but 
it has its share of complexities including the uncer- 
tainty in the initial mass function. In addition to 
the stellar mass, most galaxies have an important 
(if not dominant) fraction of their cold baryons 
in the form of neutral hydrogen gas. For consis- 
tency, in this paper we make use of stellar popula- 
tion synthesis estimates of stellar masses whenever 
possible. 

The baryonic Tully-Fisher relation (BTF) is 
one way of displaying the amount of cold baryons 
in galaxies. The BTF relation has been inves- 



tigated over the years ( McGau gh et al. _ 2000t 



Bell fc de Jo^ l200ll: l^riieiienl l200lt iMcGaugh 



20051 : IStark et all l2009t iMcGaugh et all l201 



Here we us e the r ecent observational sain ples of 
Stark et all (120091) an d lLerov et al.l (l2008l) . aloni 



with 



VerheiienI (|200ll) and the ICeha et al.l (|20& 



Stark et al 



sampl e used for the LV relation 
(|2009| ) include gas-dominated spiral galaxies, 
which makes the results much less sensitive to 
the uncertainties in the IMF. For consistency, we 



and iGeha et al.l (|2006l ) samples using a simple lin- 
ear fit to the distribution of T^-band mass-to-light 
ratios vs. (B — V) colo r shown in Figure 18 of 
Blanton fc Roweid (|2007l ): 



calculate stellar masses for the IStark et al.l ( 20091 ) 



M/Lv = 3.0(5 -V)- 0.6. (8) 

Unlike the iBell et al.l ()2003f ) models, this relation 
fits well the measured mass-to-light ratios of both 
blue and red galaxies in the SDSS. These estimates 
are fully compatible with the stellar masses used 
in the stellar mass fun ction of Li fc White! ( 20091 ) 
as part of our model. iLerov et all (|2008l) present 
results base d on the HI Nearb y Galaxy Survey 
(THINGS): I Walter et al] (|2008[ ). The measure- 
ment of luminosities in the infrared using Spitzer 
results in reliable estimates of the stellar masses 
that are consistent with the iBlanton fc Roweid 
(|2007l ) results. In this case we a djust their stellar 
masses from t he Kroupa to the Chabrieii (|2003l ) 
IMF used by iBlanton fc Roweisi (|2007l ) bv sub- 
tracting 0.05 dex. Selecting only galaxies with 
high inclinations {i > 45° or b/a > 0.7) and better 
than 15% accuracy in the circular velocity data 
leaves a total of 161 galaxies. We also include 
the result s of mass modeling of the Milky Way 
and M31 (JKlvpin et al.ll2002t IWidrow fc DubinsS 
2001 . 

We also employ equation ([8|) to obtain stellar 
masses for the early-type galaxies. This ensures 
a fair comparison between the baryonic masses of 
early- and late- type galaxies. 

Figure [TU] shows the cold baryon fraction rela- 
tive to the universal value as a function of stel- 
lar mass in our model. The cold baryon fraction 
fb peaks at w 0.2 for the stellar masses typical 
of Milky Way-type galaxies and sharply falls on 
both sides of the mass s pectrum. O u r resu lts are 
broadly consistent with IGuo et al.l (|2010l ). We 
note that even the peak of fb ~ 0.2 is almost a 
factor of two smaller than wha t a few years a go 
was considered a fiducial value (|Mo et al.lll998l ). 

The baryonic Tully-Fisher relation is shown in 
Figure [TTl Theoretical estimates from abundance 
matching provide a good fit to observational re- 
sults for galaxies ranging from dwarfs with V^irc ~ 
60 km s^^ to giants with Vdrc ~ 500 km s~^. In a 
remarkable agreement with the LV relation result, 
the model with adiabatic contraction seems to also 
provide a better fit to the BTF compared to the 
model with no contraction. There is a hint that 
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Fig. 11. — Mass in cold baryons as a function of circular velocity. The solid curve shows the median values 
for the ACDM model using halo abundance matching and including adiabatic halo contraction. The cold 
baryonic mass includes stars and cold gas and the circular velocity is measured at a galactocentric distance 
of 10 kpc. The dot-dashed curve shows the effect of neglecting halo contraction. For comparison we show 
the median and 1-ct scatter values of several binned galaxy samples. Intermediate mass galaxies such as the 
Milky Way and M31 lie very close to our model results. 



observations show more baryonic mass for dwarfs 
below Vcirc = 40 km s~^ as compared to an ex- 
trapolation of the model. It is not clear whether 
this is a real problem because of the uncertain- 
ties involved in the observations. First, the small 
sample size could produce biased results. Second, 
there is an uncertainty at the faint end of the lumi- 
nosity function. The results of abundance match- 
ing are sensitive to the number density of galax- 
ies with absolute magnitudes M^ > —14, which is 
poorly constrained. 

As in the case of the LV relation, the model 
BTF relation agrees very well with the aver- 
age population of galaxies in each morphological 



regime. Below ~ 200 km s""'^ it follows late-type 
disks while it accurately describes massive early 
types above this threshold. The observations show 
no preference for a model with no halo contrac- 
tion vs. one with maximum contraction. Both 
cases fit well within the systematic and statistical 
uncertainties in the observations. 

Although SO and elliptical galaxies seem to con- 
tain slightly less mass in cold baryons than mas- 
sive spirals, there is also a hint that the bimodal- 
ity observed in the observed LV relation in Fig- 
ures [2] and H] between early- and late-type galax- 
ies is not merely the result of a variation in the 
mass- to- light ratio. Other authors have come to 
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Fig. 10. — Baryon fraction (in stars and cold gas) 
relative to the universal value as a function of stel- 
lar mass for the ACDM model using halo abun- 
dance matching. The solid and dashed lines show 
the median and l-cr scatter of the distribution re- 
spectively. 



the same conclusio ns fe.g.. IWilliams et allbOlOt 
iDutton et al.l l2010a. bV This would imply that 
early types are not just the result of passive fading 
of late- type disks but are fundamentally different. 
It also requires that they inhabit deeper poten- 
tial wells which may be the result of different for- 
mation or environmental processes. These results 
have deep implications for galaxy formation but 
in order to draw conclusions we would need con- 
sistent stellar and gas mass estimates for a larger 
sample of galaxies, which are not currently avail- 
able. 

6.4. Galaxy Circular Velocity Function 

Projecting the distribution of galaxies in the LV 
plane onto the luminosity axis produces the lumi- 
nosity function, while projecting onto the circular 
velocity axis yields the circular velocity function 
(VF) of galaxies: the number-density of galaxies 
with given circular velocity. From a theoretical 
cosmology point of view, the VF is an ideal char- 
acterization because it does not include uncertain 
predictions for the luminosity and requires rela- 



tively modest corrections for the baryonic masses. 
Unfortunately, it is more difficult to obtain it 
from observations and so far , there have been onl> 



irom observations and so tar , tnere nave been only 
a few attempts to do s o ( Gonzalez et al.l pOOCt 



I 



Chad 



Kochanck fc White"200l HZavala et al.l l200' 
2010; Zwaan ct al. 2ol3^. 

For the theory the starting point i s the velocity 



func tion of dark matter halos (e.g., iKlvpin et al 
201Q) . For halos with Kirc < 500 km s~^ it is well 
approximated by a power-law n{> T4iic) oc V^.^", 
where a ~ 3 . This only applies to velocities taken 
at the maximum of the circular velocity curves 
of DM halos. For galaxies, the results must be 
adjusted to Vio and corrected for the dynamical 
effects of cold baryons. 

The most recent measurement of the VF 
of nearby l a te-typ e galaxies was obtained by 



Zwaan et al.l (|2010l ). Their result is based on the 
blind HI sample of the HIPASS survey, which 
is complete down to Mhi = 5.5 x 1 0^ Mq at a 
distance of 5 Mpc (|Zwaan et al.ll2010l) . Since gas- 



rich galaxies are thought to dominate at the low 
mass end, their sample should provide an accu- 
rate measurement of the abundance of dwarfs if 
these galaxies contain enough neutral gas to be de- 
tected. To obtain a galaxy velocity function for all 
morphological types we also include the determ i- 
nation of the early- type VF done bv IChad (20101), 
using the conversion betwee n velocity dispersio n 
and circular velocity found in lZwaan et al.l ( 20101 ). 
Even though their VF was obtained indirectly us- 
ing the observed relation between luminosity and 
stellar velocity dispersion, it agrees with previous 
direct measurements. 

Figure [TH shows the results, as well as the mod- 
ified Scliechte£fit_to the VF of late-type galax- 
ies ( Zwaan et al.l 20101) and the fit for early types 
found in IChad (J2010I ). At intermediate to large 
masses (80 < V^irc < 400 km s~^), where the com- 
pleteness of the surveys is hard to question, the 
VF of our model sample reproduces the observed 
abundances reasonably well. The abundance of 
MW-type galaxies is predicted to within 50% when 
adiabatic contraction is taken into account, and 
within a few percent when no contraction takes 
place. From our earlier analysis of the LV relation 
in Section IHJ we are led to believe that halo con- 
traction is needed to obtain the correct position 
of elliptical and SO galaxies in the plot. A more 
detailed treatment of AC might be necessary in 
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Fig. 12. — Comparison of theoretical (dot-dashed and thick sohd curves) and observational (dashed curve) 
circular velocity functions. The dot-dashed line shows the effect of adding the cold baryons (stellar and cold 
gas components) to the central region of each DM halo and measuring the circular velocity at 10 kpc. The 
thick solid line is the distribution obtained when the adiabatic contraction of the DM halos is considered. 
Because of uncertainties in the AC models, realistic theoretical predictions should lie between the dot- 
dashed and solid curves. Both the theory and observations are highly uncertain for rare galaxies with 
Kirc > 400 km s^^. Two vertical dotted lines divide the VF into three domains: V^hc > 400 km s~^ with 
large observational and theoretical uncertainties; < 80 km s"""^ < Vdrc < 400 km s~^ with a reasonable 
agreement, and Vi^iic < 80 km s^^, where the theory significantly overpredicts the abundance of dwarfs. 



order to better match the abundance of galaxies 
larger than the Milky Way. Our model galaxy 
VF overestimates the abundance of the most mas- 
sive and rarest galaxies with VcItc > 400 km s^^ 
regardless of whether or not we implement the 
correction for contraction of the halos. Most of 
these extremely bright galaxies inhabit the cen- 
ters of clusters, where it is very likely that the 
simplistic observational estimate of l^ii-c is break- 
ing down. At small velocities {Vdrc < 80 km s^^) 
the theory significantly overpredicts the number of 



dwarfs. This "missing dwarfs" problem remains 
unresol ved in ACDM dTikhonoy fc KlvpinI l2009t 



Zavala et all 120091 : IZwaan et al.ll2010[ ). It should 



be noted that the variance of the velocity function 
of the model galaxies below 60 km s~^ in regions 
of radius 5 Mpc can be as large as 1 order of mag- 
nitude. This shows that environmental bias may 
be an important factor in explaining the under- 
abundance of dwarfs in our model compared to 
HIPASS. 

To illustrate the effect that each of the steps 
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Fig. 13. — Details of the velocity function. The 
dotted line corresponds to the dark matter halo 
VF at z = 0, while the thick solid line shows the 
distribution of galaxies obtained if the maximum 
rotation velocity of the halos is measured at its 
historical maximum (i.e. before accretion). Note 
that the total mass includes 17% in baryons that 
behave like dark matter in dissipationless simula- 
tions. The short (long) dashed curve shows the 
Schechter fit for late (early) type galaxies. The 
thin full curve is the total observed VF. 

in our procedure has on the VF, we show in Fig- 
ure [13] the VF of DM halos only. It also shows 
that when the stripping due to the merger his- 
tory of each halo is considered, the halo VF does 
a slightly better job at matching the abundance of 
galaxies. 

As we previously noted, the corrections due to 
the presence of the cold baryonic component affect 
dwarfs (Vchc < 100 km s~-^) very little, resulting 
in a negligible shift in their abundance compared 
to that of their host DM halos at the low-mass 
end of the VF in Figure [T^ One interpretation 
of this is that the dwarf overabundance problem 
cannot be resolved if both the LV relation and the 
VF of dwarf galaxies are to be reproduced simul- 
taneously. In other words, these galaxies must un- 
dergo a process that limits their abundance with- 
out changing their dynamical mass. The first pos- 
sible origin for the large discrepancy between our 



model galaxies and the HIPASS VF could be ob- 
servational bias. HIPASS is a blind HI survey and 
does not detect gas-poor galaxies. Only if gas- 
poor dwarf spheroidals dominate the galaxy pop- 
ulation below ~ 100 km s~^ would it be possible 
to reconcile our results with the survey. This is 
highly unlikely since this type of galaxies are only 
a small fraction of the total dwarf population. On 
the other hand, if the HIPASS HI mass detection 
limit (5.5 x lO'' M© at 5 Mpc) is relatively high at 
the distances where most of their sample is found, 
incompleteness effects might explain the discrep- 
ancy. 

Assuming that the surveys are complete, a pos- 
sible solution to the problem is a mapping of all 
the dwarf galaxies below 50 km s~^ to DM halos 
in the range 50 — 100 km s~^. This in turn im- 
plies that the measured rotation curves of a large 
fraction of dwarfs must severely underestimate the 
true maximum circular velocities of these galaxies. 
The only possible explanation for this bias would 
be that the optical and HI disk is truncated well 
inside the radius where the rotation curve flattens 
out. Another solution to the missing dwarf prob- 
lem requires most of these galaxies to have a low 
enough surface brightness in HI to be undetectable 
in current surveys. This would imply the existence 
of a large number of small halos containing little 
or no neutral gas. 

6.5. Galaxy two-point correlation function 

The most important success of the halo abun- 
dance matching technique is considered to be re- 
producing the observed galaxy clustering mea- 
sured in the form of the galaxy correlation func- 
tion in its various forms, both in the local uni- 
verse and at high rcdshift (Tasitsiomi ct al. 200 



Conrov et al.lbOOfiUCuo et aLbOldlWetzel fc Whit: 



20101 ). Most of these works claimed to match the 
observed clustering although they relied on simu- 
lations with either very low resolution or outdated 
cosmological parameters. The high resolution and 
large volume of the Bolshoi simulation allow us 
to calculate the galaxy two-point correlation func- 
tion at a range of scales comparable to the latest 
results from the fin al data release of the SDSS 
(jZehavi et al.ll2010l) . In addition, its up-to-date 
set of cosmological parameters allows for a direct 
comparison between observations and the predic- 
tions of ACDM+HAM. The comparisons in this 
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section do not make use of the dynamical correc- 
tions that were necessary to obtain the LV relation 
and the velocity function. Instead, the calculation 
of the galaxy correlation function only requires 
the position and velocity information of the ha- 
los in the simulation along with their luminosities 
obtained from HAM. This makes the correlation 
function an even more robust prediction of the 
model. 

In order to compare our model with observa- 
tions, we use the most recent measurement of the 
SDSS g alaxy projected aut ocorrelation function 
done by IZehavi et al.l ( 20101 ). To make the best 
comparison possible we use projected galaxy sep- 
arations (a direct observable) and the same lu- 
minos ity and projected radii bins as IZehavi et al 



( 2010( ). We also integrate along the line of sight 
using the same distance bins while including the 
peculiar velocities of the model galaxies in the red- 
shift calculation. The integration is traditionally 
performed to wash out the effects of redshift dis- 
tortions. We limit the calculation of the corre- 
lation function to distances below 30 h^^ Mpc to 
avoid scales at which much of the power comes 
from long waves that are absent in the simulation 
due to the finite box size. The small-scale correla- 
tion function of DM halos is extremely sensitive to 
the abundance of satellite halos near the centers 
of hosts. As a result of this, the clustering in A^- 
body simulations could be underestimated due to 
artificial disruption of just a few satellites. Since 
it is beyond the scope of this paper to perform 
a comprehensive study of this effect, we choose 
to compare our model to galaxies in the range 
-19 > Mr - 5 log /i > -22. In this and aU follow- 
ing sections we refer to the "'^r-band in shorthand 
as simply the r-band. 

Figure [14] shows the projected two-point cor- 
relation function of the model galaxies in Bolshoi 
and compares it to the full SDSS sample results. 
The clustering amplitude of the model galaxies 
is in excellent agreement with observations for 
galaxies with luminosities around L* in the range 
-20 > Mr - 5 logio h> -21. Model galaxies with 
luminosities —19 > Mr — Slogj^g''- > ~21 agree 
very well with the observations at scales beyond 
1 — 2 h"^ Mpc where the clustering is dominated 
by halos of different hosts (the so-called two-halo 
term). Below 1 h^^ Mpc there is a marked decline 
in the number of pairs as the separation decreases. 



For bright galaxies with — 21 > Mr — Slogj^p ^^ > 
-22 the situation is different; ACDM + HAM 
slightly overpredicts the clustering over all scales, 
with the disagreement increasing to ^ 30% beyond 
10 h-^ Mpc. 

The discrepancy in the clustering of the faintest 
bin at small separations may be a result of nu- 
merical effects such as artificial disruption or 
halo misidentification in dense environments. 
A small deviation at the closest separations 
{rp < 400 h"^ kpc) is likely to have the same 
origin. Since the fraction of halos that are satel- 
lites decreases sha rply at large halo masses (see 
Klvpin et al.l 120101 1 , the model correlation func- 
tions of the brightest galaxies do not suffer from 
these effects. Further scrutiny is necessary to un- 
derstand the origin of the effect and make more 
robust comparisons with observations. 

6.5.1. Effect of scatter on the correlation func- 
tion 

Since only the brightest galaxies in the LV rela- 
tion are affected by scatter, the obtained two- 
point correlation function of galaxies brighter 
than Mr ~ —22 will be sensitive to the choice 
of scatter distribution. In the past few years, 
some studies of the correlation function of DM 
halos have suggested the scatter to be an es- 
sential in gredient in reproducing the observa- 
tions (e.g..lTasitsiomi eraLl l2004l: I Wetzel fc Whit3 



2OIOI : iBehroozi et al.ll2oTol ). 



Figure [15] shows the galaxy autocorrelation 
function obtained for our model galaxies using 
the stochastic method described in Appendix [A] 
to perform HAM. Here we assume the same distri- 
bution described in Section [6.2.41 with cta/^ « 0.5 
below L* and ctm^ ~ 0.3 above. The correla- 
tion amplitude of model galaxies fainter than 
Mr— 5 logj^o /i = — 21 is essentially unchanged com- 
pared to the monotonic assignment result shown 
in Figure [Ml The clustering of the bright galaxies 
with —22 < Mr — 51og]^o ^ < ~21 shows a slight 
decrease at all separations except the smallest 
ones, and thus better agreement with the SDSS 
data. The amplitude decreased due to the fact 
that the same galaxies now get assigned to less 
massive halos on average, and these halos are less 
clustered. This is consistent with the upward shift 
in the average luminosity of the brightest galax- 
ies in the LV relation (Figure [9]). The correlation 
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Fig. 14. — The correlation function of the Bolshoi galaxies using HAM without scatter vs. the SDSS 
observations. Top left: correlation function of model galaxies in three magnitude bins showing the poisson 
uncertainties as thin lines. Top right, Bottom left, Bottom right: the clus tering in each lumin osity bin is 
compared to SDSS galaxies. Solid circles with error bars are the data from lZehavi et al.l (|201Cl( ). ACDM + 
HAM does an excellent job at reproducing the shape and amplitude of the clustering of galaxies near the 
the knee of the luminosity function (—20 > Mr — 5 logj^p ^ > ~21). Brighter model galaxies are slightly more 
clustered than SDSS galaxies at large separations while faint ones underestimate the observed clustering at 
distances below 0.5 h~^ Mpc. 



functions of galaxies in fainter bins are indistin- 
guishable from those without scatter. 

In summary, applying our physically motivated 



scatter model maintains (and even improves) the 
excellent agreement of the ACDM+HAM model 
with the observed galaxy clustering. The cluster- 
ing of the most massive and rare galaxies (those 
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Fig. 15. — Same as Figure [Ml but including variable scatter in luminosity at fixed circular velocity using our 
stochastic abundance matching method. Galaxies with —19 > M,. — Slogj^p/i > —21 are mostly unaffected 
while those in the brightest bin are slightly less clustered at all separations than in the case with no scatter. 
This results in a better agreement with the SDSS observations. 



above Air ~ —22) will be more sensitive to the ad- 
dition of scatter and the model used to implement 
it. Once other dominant sources of uncertainty 
in the simulations and observations are better un- 
derstood, a robust test of the cosmological model 
could be done using these objects. 



7. Comparison v^rith other results 



Ou r results are broadly consistent with lGuo et al 
( 20101 ). who also use the abundance matching 
technique. Specifically, in their Figure 6 they 
show the stellar-mass - circular velocity rela- 
tion. The theoretical velocities appear to be 
smaller than the observed circular velocities for 
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100 - 150 km s~^ Although iGuo et al 



(J2OIOI ) did not apply the necessary corrections 
discussed in our paper, they argue that inclusion 
of the cold baryon mass may bring the theory into 
agreement with the observations. As we show, 
indeed this is the case. 

I ncidentally, in the semi-analytic modeling pa- 
per (JGuo et al.l 120111) based on the MiUennium-I 
and II simulations, the predicted angular correla- 
tion function of galaxies with logAf* < 10.77 is 
significantly too high compared with SDSS data, 
especially at separations less than about 1 Mpc. 
The authors attribute this to the fact that the 
large CTg = 0.90 used in the Millennium simula- 
tions produced too many massive halos that in 
turn host too many pairs of galaxies in their sub- 
halos. 



Button et al.l (|2007l ) argue that the standard 



cosmological model with adiabatic contraction and 
standard concentrations fails to simultaneously 
reproduce the observed LV relation and the lu- 
minosity function for late-type galaxies. This 
conclusion is not compatible with our results. 



A num ber of assumptions made in iDutton et al 
(|2007l ) are either outdated or need corrections. 



For example, for their preferred m odel they use 
the "standard concentrations" of Bullock et al. 



(J2001I ). which were based on a simulation with 
(Tg = 0.9, although they attempted to rescale them 
to a cosmological model with the normalization 
a% = 0.8. The normalization of the current ACDM 
cosmological model is as = 0.82 base d on CMB 
and other data (e.g. ljarosik et al.ll2010l ). which re- 
sults in halo concentrations that are ~30% lower 



than what iDutton et al.l ([20071) used in their pre- 
ferred model. In turn, this reduces the dark mat- 
ter circular velocities in the inner regions of ha- 
los by about 15%. Som e of the neces s ary co rrec- 
tions were discussed bv lOutton et al.l (|2007l ) and 
it was shown that they substantially improve the 
fit of the TF relation. However, the main dif- 
ferenc e is the treatment o f the luminosity func- 
tion. iDutton et al.l (|2007l ) use criteria obtained 
from SAMs to argue that a model with halo ex- 
pansion is necessary to match the LF. We avoid 
such assumptions completely because our model 
reproduces the galaxy statistics automatically. 



Gnedin et al.l (|2007D studied structural proper- 



(|200l[ ) high concentrations as the "standard" 
model. It was concluded that the theory has 
problems and that adiabatic contraction is the 
likely culprit. Another possible solution was to 
lower the ha l o coii centrations. Indeed, when 
iGnedin et al.l (J2007l) used concentrations for a 
model with a% = 0.74 as predicted by simple the- 
oretical arguments, they found that the theory 
gives an acceptable fit to the data. The prob- 
lem is that a^ = 0.74 is too low. However, it 
seems that their analytical scaling with ag, was 
not accurate enough: the concentrations actu- 
ally used by iGnedin et al.l ([20071) are practically 
(within 3%) the same as what we find in iV-body 
simulations for the Bolshoi A CDM model with 
as. = 0.82 (JKlvpin et al.ll2010t ). In short, there 
seems to be no contradiction between our results 



ties of spiral galaxies and compared t hem with the- 
oretical predictions. They also used lBuUock et al 



and IGnedin et al.l (|2007l ) even when we consider 
models with standard adiabatic contraction. More 
definite conclusions require careful analysis and 
changes in the fraction of cold baryons among 
other things. 

8. Discussion 

In this paper we address one of the most diffi- 
cult problems in cosmology: is the standard cos- 
mological ACDM model compatible with observa- 
tions when it comes to the prediction of the abun- 
dance and properties of galaxies? Instead of fo- 
cusing on traditional issues such as the zero-point 
and the slope of the Tully-Fisher relation for spiral 
galaxies, we work with a more generic luminosity- 
velocity (LV) relation: a correlation of galaxy lu- 
minosity with the circular velocity at a 10 kpc ra- 
dius. We also investigate the (cold) baryonic mass 
- velocity relation, which following tradition we 
call the baryonic Tully-Fisher (BTF) relation, as 
well as the velocity function and the two-point au- 
tocorrelation function of galaxies. All these statis- 
tics encompass galaxies of different types - from 
dwarf galaxies to normal spirals to giant ellipticals. 
These statistics - in combination with the theo- 
retical predictions of the cosmic microwave back- 
ground and the abundance and properties of dark 
matter halos - are major tests for the validity of 
the ACDM model. 

We use the abundance matching technique to 
assign luminosities to halos predicted by cosmolog- 
ical simulations. We also use abundance matching 
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to assign stellar and cold baryon masses. We find 
that all three statistics - the LV and BTF rela- 
tions, and the velocity function - provide reason- 
ably good fits to observations for galaxies ranging 
over 10 magnitudes in luminosity and for circular 
velocities from 80 km s^^ to 400 km s~^. By con- 
struction, our models fit the observed luminosity 
and stellar mass functio ns. Since they are ba sed on 
the Bolshoi simulation (jKlvpin et al.l 120101 ) . they 
also fit known properties of dark matter halos in- 
cluding the halo mass function and the dependence 
of halo concentration on mass. In addition, in this 
paper we show that halo abundance matching also 
yields the correct clustering properties of bright 
galaxies. In short, we have a model, that fits - at 
least on average - all the basic statistics of galax- 
ies with T4irc > 80 km s~^ considered at a ^10 kpc 
scale. 

Matching theory with observations requires a 
careful consideration of many different effects and 
application of different corrections. These effects 
were considered both for observations and for the 
theory. On the observational side, we compiled a 
representative sample of galaxies with measured 
circular velocities. Velocities were either asymp- 
totic values ("flat part" of rotation curves) for 
spirals or measurements at ^ 10 kpc radius for 
SOs and Es. We do not use fits (such as power- 
laws) to the data but instead work directly with 
the distributions. We do not apply morphological 
corrections of the TF relation (e.g., differences be- 
tween Sb and Sa galaxies) because those corrupt 
the bright end of the LV relation. Since the Tully- 
Fisher luminosities are corrected to face-on, we de- 
correct the magnitudes of galaxies for the effect of 
internal absorption to make them consistent with 
the measurement of the luminosity function. 

For the theoretical predictions we try to make 
all the possible corrections to mimic the observa- 
tional situation. For example, we do not use virial 
masses of halos because virial radii are too large 
compared with the typical distances at which ro- 
tational velocities of observed galaxies are mea- 
sured. We do not assume a particular shape for 
the halo density profiles: they are measured di- 
rectly in the simulations. The simulations required 
for this type of analysis should have a very high 
resolution so that subhalos are also resolved. This 
allows us to avoid using intermediate steps such 
as the Halo Occupation Distribution or the Con- 



ditional Luminosity Function, which are often ap- 
plied to low-resolution simul ations. T he Bolshoi 
cosmological simulation (Klv pin et al.l 12010) pro- 
vides high quality results resolving distinct halos 
and subhalos down to the completeness limit of 
T4irc = 50 km s^^. 

The observations should be taken cautiously 
since each has a different degree of accuracy. The 
LV relation is the most accurate because it is eas- 
ier to measure luminosities than to estimate stel- 
lar masses, which require additional modeling and 
assumptions. This is why we consider the LV re- 
lation as our prime target. The velocity function 
is the least reliable since observations are still at 
the very early stages. The completeness of the 
HIPASS VF is very uncertain. Just the fact that 
the detection limit is quoted at 5 Mpc shows that 
the accuracy of the HI mass function is not very 
high. This is why we treat the results on the ve- 
locity function for V^circ > 80 km s~^ as a "pass" 
for the theory in spite of some deviations such as 
at 130 km s^^. More accurate treatment of these 
gas-rich galaxies may also change the situation: 
after all, changes in abundances and velocities by 
^ 10% may (or may not) resolve the discrepancies. 

It is more difficult to reconcile the theory and 
observations at smaller velocities. Indeed, at 
V^irc — 50 km s^^ the formal disagreement is al- 
most a factor of ten. This is the only serious prob- 
lem that we find when matching galaxies with dark 
matter halos. A similar pr oblem on somewhat 
smalle r scales was reported bv lTikhonov fc Klvpin 
(J2009I ). who studied the population of dwarfs in 
the ~ 1 Mpc region centered on the Milky Way 
galaxy. iTikhonov et al.l (|2009) argue that Warm 
Dark matter may be the solution to the problem. 

We introduce a simple scatter model that is well 
motivated and preserves the agreement with the 
LV relation and the correlation function of galax- 
ies in the SDSS. The introduction of scatter has 
some complications. Observed deviations from the 
median relations seem to have a systematic com- 
ponent: early- type galaxies are systematically be- 
low the median LV relation and gas-rich spirals 
are above it. It seems likely that the LV relation - 
like the color-magnitude diagram - has a bimodal 
structure. In this case, no simple gaussian spread 
can explain the whole diagram. However uncer- 
tain, the spread must be explained. One approach 
might be to match halos separately to red and blue 
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galaxies, for example using local density as well as 
luminosity. Ultimately it will be necessary to find 
the real source of the dynamical bimodality and 
to measure it observationally. 

Although our ACDM+HAM prescription makes 
simplifying assumptions regarding the distribution 
of baryons in DM halos, it yields results that are 
compatible with more detailed dynamical mod- 
els. Our model pr edicts Vdrc values that differ by 
less than 5% from iDutton et al.l (J2010bl ) for mas- 
sive dis k galaxies w i thout h alo contraction. Given 
that iDutton et al.l (|2010bf ) include a large set of 
observational constraints on the radial distribu- 
tion of baryons, the consistency with our results 
is evidence of the robustness of our modeO. 

Abundance matching is a very successful way 
to make predictions about how on average galax- 
ies can inhabit dark matter halos. It gives up 
solving the most difficult and the most important 
problem: how galaxies form inside dark matter ha- 
los. It simply assumes that the stellar mass and 
luminosity monotonically (or possibly with some 
scatter) scale with the circular velocity. Bluntly 
speaking, it assumes that the maximum circular 
velocity of a halo determines the properties of the 
galaxy hosted by that halo. Remarkably, this can 
reproduce some basic environmental relations such 
as the morphology-density relation and the depen- 
dence of gala xy clustering on th e luminosity of 
galaxies fe.g.. IConrov et al.ll2006( ) because of the 
correl ation of environment with the average hal o 



mass ( Sheth et al.l 120011: ISheth fc TormenI 12004 ) 



However, there are potential issues with abun- 
dance matching. It is not clear how it can explain 
dependencies on environment even if galaxies are 
selected with the s ame r-band luminosity or the 
same stellar mass (jHogg et al.l 120041 : Ivan der Wei 
20081 ). Modehng of the bimodality in the LV re- 



lation (the apparent differences between early and 
late type galaxies) is another problem to address. 



iDutton et al.l ll2010bl ) used LV relations to find that late- 
and early-type galaxies are best fitted by h alo expansion 
and halo contraction (a la lGnedin et al.l2004) , respectively. 
As mentioned above, we do see a systematic difference be- 
tween early- and late-types in our assembled data sets that 
might be explained by some combination of differences in 
halo contraction and halo masses (Figures [4l and [TTIl . How- 
ever, given the observational uncertainties and the non- 
differentiation between early and late- types in our models, 
we cannot yet provide a precise interpretation of the bi- 
modality. 



It will be interesting to see how much better the 
results will be from more sophisticated abundance 
matching including g alaxy color and local density 
- as, e.g., in ITasitsiomi et al.l ( 20041 ) - and from 
semi-analytic modeling based on the Bolshoi sim- 
ulation. This work is in progress. 

Disk formation and semi-analytic models still 
struggle to simultaneously reproduce the TF re- 
lation (a subset of the LV relation) and the 



abundance of galax i es fe.g.. iBenson et al.M2003 ; 



Monaco et all l2007t iBenson &: Boweij l2010t ). In 



this paper we have shown that our model is suc- 
cessful at this task. Simultaneously reproducing 
the luminosity function and the LV relation de- 
pends critically on implementing each of the steps 
in Section[5]to obtain the properties of the galaxies 
that inhabit ACDM halos. For example. Figure |H] 
shows how assuming an incorrect val ue for the 
baryon fraction (as in IMo et al.l ( 19981) ) can lead 
to an LV relation that is in striking disagreement 
with the observa tions. In the more rece nt semi- 
analytic model of IBenson fc BoweiJ ( 20101) the cir- 
cular velocities of galaxies are 40 — 50 km s~^ 
larger than observed at any luminosity. This may 
be explained by the fact that their baryon fraction 
is about 20% larger than ours for Milky Way-mass 
galaxies and about an order of magnitude larger 
than our result for the most massive ellipticals as 
well as dwarfs. Our model shows that galaxies 
with masses larger that the Milky Way have cir- 
cular velocities that are extremely sensitive to the 
baryon content within their optical radius which 
may explain why SAMs overpredict the circular 
velocity. 

Previous works based on abundance match- 
ing were successful at reproducing the statistics 
of the integrated properties of galaxies (such as 
clustering as a function of luminosity and red- 
shift) but made no attempt to incl ude their inter 
nal and baryonic structu re 
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Kravtsov et al 



Behroozi et al...2010.: .Guo et al...2010[ ). In this pa 



2004 IConro v et al.' '2006: 'Wetzel fc Whitd l2010t 



per we show that making robust dynamical cor- 
rections to the structure of halos obtained in sim- 
ulations gives the correct galaxy scaling relations. 
These corrections include adding the cold baryonic 
component and measuring ycirc at 10 kpc. Previ- 
ous studies using HAM could not include these 
corrections partly because simulations lacked the 
large dynamic range necessary to form the largest 
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halos and resolve substructure adequately. The 
large box size and very high resolution of the Bol- 
shoi simulation makes it possible to obtain good 
statistics of even the largest clusters and resolve 
the structure of dwarf halos. 

9. Conclusions 

Here is a short summary of our results: 

• In combination with previous results, we 
conclude that the standard ACDM model in 
conjunction with halo abundance matching 
can simultaneously fit reasonably well the 
main global statistics of galaxies: the lu- 
minosity function, the stellar mass function, 
the Luminosity- Velocity relation, the bary- 
onic Tully-Fisher relation, the abundance 
of galaxies with circular velocities T4irc > 
80 km s~^, as well as the clustering prop- 
erties of bright galaxies as a function of lu- 
minosity. 

• There are systematic deviations in the LV re- 
lation with SO and elliptical galaxies located 
about 1 magnitude below late types in the 
LV relation. Massive early types contain less 
baryonic mass than late types at the same 
circular velocity. 

• The range of the effect of contraction of the 
DM halos due to baryon infall brackets the 
observations. The LV relation shows prefer- 
ence for a model with moderate contraction, 



as predicted bv iGnedin et al.l (|2004[ ). This 
result is compatible with the observed ve- 
locity function of galaxies. 

• There seems to be an overabundance of 
model galaxies by a factor of ~ 10 compared 
to observed dwarf galaxies with Vcirc < 
50 km s~^. This is a serious problem for the 
ACDM model: galaxies with these circular 
velocities cannot be affected much by "nor- 
mal" physical processes (e.g., supernovae 
feedback or reionization of the Universe) 
proposed for the solution of the satellite 
problem at Vdrc ^ 30 km s~^. However, 
the observational results on the abundance 
of dwarf galaxies still need to be improved. 

• Including scatter in luminosity at fixed Vchc 
using our physically motivated scatter model 



maintains the agreement of the model LV 
relation with observations. 

The correlation function of the model Bol- 
shoi galaxies matches very well the observa- 
tions of bright galaxies. The agreement im- 
proves when implementing our scatter pre- 
scription for all but the brightest galaxies, 
where a better understanding of the uncer- 
tainties is necessary to make a fair test of 
ACDM using halo abundance matching. 
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APPENDIX 



A. Halo abundance matching including stochastic scatter 

To perform halo abundance matching while including scatter in luminosity as a function of circular velocity 
we perform the following procedure. For brevity, wherever we use Mr we refer to Mo.ij. — Slogj^Q h. 

1. Start with the monotonic assignment described in Section [5j Order the list of DM halos from largest 
to smallest Khc- Using the integral LF, solve for the unique luminosities of the galaxies that have 
the same number densities as the DM halos in Bolshoi. This matching gives the monotonic relation 

M-^-Cl/eirc). 

2. For the halo with the largest Vcirc, draw a luminosity value M^ at random from a Gaussian distribution 
of width a centered at a point l-cr brighter than the value of M™°"° for that halo. Mathematically, this 
is equivalent to Mr — M^°'^° {Vdrc) — <J + ^(0, a), where Q{0, a) is a random reahzation of a Gaussian 
probability distribution with standard deviation a centered at zero. 

3. If the randomly drawn Mr is brighter than the brightest galaxy in the LF, another draw is performed 
and the process is repeated until a suitable value is found. 

4. The random Mr is compared to the list of available Mr values in the list. The closest available value 
becomes the luminosity of that halo. If the value is already taken, another random draw is performed 
until an available one is found. This luminosity value is flagged to prevent it from being used again for 
another halo. This step ensures that the observed luminosity function is preserved by only assigning 
each luminosity once. 

5. In order to avoid having unassigned values of Mr in the list on the bright tail of the distribution, we 
check whether there is any unassigned luminosity which is more than 3-cr brighter than M™°"°. If such 
unassigned value exists, the next halo in the list is assigned to it. Since we are stepping along the list 
of luminosities as we assign them, the process is intrinsically asymmetric and there will always be some 
leftover Mr values that get assigned in this step. The offset used in step 1 ensures that these make up 
only a few percent of the sample. 

6. Repeat steps 2-5 for each halo in the ordered list until the 3-cr faint tail of the Gaussian for a given 
halo reaches the completeness limit of the sample as defined in Section [3l The procedure is stopped at 
this point to prevent placing a hard constraint on the faint end of the LF where it is most uncertain. 
The faintest halos (up to 0.5 magnitudes brighter than the cutoff) are removed from the sample to 
insure the preservation of the LF throughout. 

In spite of the fact that a constant Gaussian distribution is used in the algorithm, the final distribution 
of Mr for a given circular velocity Vi^irc is not a Gaussian and the width of the Mr — Vi^irc relation is not 
constant. This happens because of the asymmetry in the distribution of the galaxies: there are always more 
galaxies with smaller luminosities than with larger ones. Thus, the distribution and the spread of the final 
Mr — Vcirc relation are the result of a convolution of a Gaussian distribution with the luminosity function. 

The effect of the asymmetry of the LF is more pronounced for the brightest galaxies in the sample and 
accounts for the small shift in the median compared to the monotonic assignment. The algorithm also 
gives a natural narrowing of the distribution of luminosities as T^irc increases. This occurs because in the 
exponential tail of the LF the number of available Mr values changes rapidly across the width of the Gaussian. 
Increasingly fewer available values on the bright side of the median force the selection of most values to take 
place in a narrower interval on the faint side. 
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Since our assignment method reduces the width of the obtained distribution of luminosities, we choose the 
input value a = 0.7. This yields a distribution with gm^ ~ 0.5 below ~ 200 km s~^ and gradually decreasing 
to (JM^ ~ 0.3 above ~ 300 km s^^. 

Figure [TCI shows the distribution of luminosities in four bins of circular velocity, Viq. Below 200 km s~^, 
galaxies show a near-normal distribution that is centered very close to the values obtained from the monotonic 
assignment without scatter. For 14irc > 250 km s~^ the distributions get progressively more skewed as 
galaxies move from the bright to the faint tail. The spread of the distribution of luminosities increases 
with V^ii-c^ the widths are aMr ~ 0.50, 0.45, 0.43 and 0.35 for the bins centered at 102.5, 205.0, 307.5 and 
520.0 km s^^ respectively. 







400 



300 



200 



100 - 



-20 -21 -22 -23 




21 -22 







2000 



1000 - 




19 -20 -21 



■16 -17 -18 -19 



Mo,^ - Slog^oh 



Fig. 16. — The distribution of model galaxies obtained using the stochastic abundance matching method. 
Each panel shows one of four representative circular velocity (Vio) bins. The vertical dashed lines in each 
panel show the median and average while the vertical solid line shows the average value that was assigned in 
the monotonic scheme. The dotted lines show gaussian fits to each distribution. As galaxies become brighter 
(from bottom right to top left), the distribution narrows and becomes slightly skewed towards the faint tail. 



Figure [IT] shows the LV relation obtained with the scatter model of iBehroozi et al.l (20101) in the case of 
constant scatter (left panel) as well as assuming the same variable width used in our model: gm^ ~ 0.5 for 
14irc < 250 km s^^ declining to ^ 0.3 for larger Kirc- Evidently, the median of the distribution in luminosity 
stays relatively unchanged only when forcing a small scatter width at the bright end. Even when we impose 
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the same variable width used in the stochastic HAM scheme, the resuhing spread is too large for the brightest 
galaxies. 
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Fig. 17. — The LV relation of the Bolshoi galaxies obtained using the deconvolution method assuming log- 
normal scatter. The solid (dotted) line shows the median (average) of the circular velocity in bins of r-band 
luminosity. The shaded areas encompass 68% and 95% of the galaxies in each bin. The dashed line shows 
the result of monotonic assignment with no scatter. Left: result of using a constant scatter width aMr = 0.5. 
The median relation deviates by up to 1 magnitude for the brightest galaxies compared with the monotonic 
result shown as a dashed line. Right: result of using a width aM^ = 0.5 below L* and ctm,. = 0.3 above. 
Although the median only deviates by a small amount, the spread in luminosity of the distribution of galaxies 
with Mr — 51og]^Q h < —20 is considerably larger than that obtained using stochastic HAM. 



B. Early- type data 

Table [B] gives some of the properties of our early-type sample along with the source of data for each 
galaxy. 
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Table 1 
Luminosity and circular velocity data for nearby early- type galaxies. 



Name 



Type 



Mb 



l°Sio {^] 



Vio 



+ 



Vl5 



+ 



V20 



+ 



probe 



ref. 



;C 4889 


E 


-22.61 


0.16 


11.87 


520 


60 


;C 4874 


E 


-22.51 


0.16 


11.65 


500 


30 


;C 0315 


E 


-22.27 


0.25 


11.66 


520 


10 


;C 1316 


SO 


-22.20 


0.19 


11.46 


381 


29 


;C 4839 


E 


-22.14 


0.15 


11.52 


385 


65 


;C 0057 


E 


-22.06 


0.19 


11.56 


491 


20 


;C 4555 


E 


-22.05 


0.21 


11.50 


614 


58 


;C 4952 


E 


-21.75 


0.19 


11.37 


405 


10 


;C 6407 


E/SO 


-21.73 


0.26 


11.51 


445 


25 


;C 4472 


E 


-21.64 


0.12 


11.41 


415 


40 


;C 7626 


E 


-21.56 


0.20 


11.41 


420 


9 


;C 5044 


E 


-21.49 


0.27 


11.31 


273 


17 


;C 4486 


E 


-21.40 


0.16 


11.31 


503 


47 


;C 3923 


E 


-21.35 


0.43 


11.26 


365 


28 


;C 4816 


E/SO 


-21.35 


0.23 


11.19 


300 


45 


;C 1395 


E 


-21.32 


0.17 


11.26 


374 


17 


;C 4944 


SO 


-21.31 


0.32 


11.19 


275 


2 


;C 4382 


SO/a 


-21.31 


0.16 


11.12 


260 


19 


;C 4649 


E 


-21.29 


0.16 


11.27 


425 


10 


;C 4374 


E 


-21.25 


0.12 


11.26 


382 


20 


;C 4827 


E/SO 


-21.25 


0.23 


11.20 


350 


50 


;C 0128 


SO 


-21.20 


0.20 


11.18 


370 


14 


1459 


E 


-21.17 


0.21 


11.26 


338 


53 


;C 4957 


E 


-21.14 


0.19 


11.21 


325 


2 


;C 4261 


E 


-21.04 


0.20 


11.20 


362 


27 


;C 6703 


E/SO 


-21.07 


0.20 


11.11 


220 


20 


;C 3665 


SO 


-21.06 


0.23 


11.13 


423 


52 


;C 5846 


E 


-21.04 


0.24 


11.19 


340 


5 


;C 4908 


E 


-21.00 


0.15 


11.19 


320 


40 


;C 0720 


E 


-20.99 


0.18 


11.18 


317 


12 


;C 7796 


E 


-20.99 


0.20 


11.15 


308 


46 


;C 4365 


E 


-20.98 


0.18 


11.15 


333 


26 


;C 3607 


E/SO 


-20.91 


0.20 


11.07 


278 


28 


;C 1399 


E 


-20.88 


0.19 


11.11 


430 


25 


;C 3585 


E 


-20.77 


0.22 


10.99 


295 


46 


;C 2974 


E 


-20.76 


0.20 


11.06 


304 


10 


;C 5084 


SO 


-20.73 


0.21 


11.05 


282 


8 


;C 6771 


SO/a 


-20.72 


0.17 


10.92 


340 


18 


;C 4807 


E/SO 


-20.64 


0.23 


10.99 


295 


30 


;C 4931 


SO 


-20.62 


0.23 


11.02 


280 


5 


;C 0821 


E 


-20.58 


0.21 


10.88 


182 


13 


;C 1332 


E/SO 


-20.56 


0.22 


10.94 


291 


9 


0843 


SO 


-20.55 


0.19 


10.83 


380 


10 


O151-G004 


SO 


-20.48 


0.26 


11.10 


291 


19 


;C 4494 


E 


-20.40 


0.17 


10.77 


198 


10 


;C 4869 


E 


-20.38 


0.15 


10.91 


280 


50 


;C 4636 


E 


-20.38 


0.17 


10.85 


430 


16 


;C 1032 


SO/a 


-20.33 


0.21 


10.82 


270 


20 


;C 4697 


E 


-20.20 


0.18 


10.73 


235 


15 


4045 


E 


-20.19 


0.15 


10.85 


390 


40 



50 — — — — — — stars T-l-07 

50 510 60 70 — — — stars T-l-07 

5 530 30 10 550 35 20 stars K-fOO 

29 466 52 52 361 49 49 X-ray NM09 

30 400 80 50 425 105 70 stars T-l-07 
49 494 19 52 495 22 52 X-ray O-f-07 
62 606 61 66 598 63 74 X-ray OP04 
15 435 20 30 455 25 45 stars T-l-07 
35 450 50 40 460 60 60 stars MBOl 
40 — — — — — — stars MBOl 

10 410 15 15 390 20 15 stars K-I-OO 

17 328 11 11 365 13 13 X-ray NM09 

34 — — — — — — GCs M-l-11 

28 324 40 40 324 40 40 X-ray NM09 
40 300 50 40 310 50 40 stars T-l-07 

17 — — — — — — X-ray NM09 

2 280 2 2 280 2 2 stars T+07 

19 260 19 19 198 51 51 X-ray NM09 
10 436 21 22 444 35 37 X-ray HBIO 

20 386 20 20 393 20 20 PNe N-l-11 
30 — — — — — — stars T-l-07 

14 361 15 15 — — — stars W-f09 
53 258 53 53 282 45 45 X-ray NM09 
2 310 10 2 295 15 2 stars T-l-07 

29 338 27 29 332 30 33 X-ray HBIO 

15 — — — — — — stars K-I-OO 

52 404 32 32 404 32 32 X-ray NM09 

5 — — — — — — stars K-fOO 

40 — — — — — — stars T+07 

13 323 13 13 330 11 12 X-ray HBIO 

26 305 37 28 297 35 26 X-ray O+07 

26 333 26 26 — — — X-ray NM09 

28 278 28 28 265 28 28 X-ray NM09 

30 — — — — — — stars K-I-OO 

46 — — — — — — X-ray NM09 

10 — — — — — — gas W-l-08 

8 — — — — — — stars W-l-09 

18 331 18 18 320 23 23 stars W+09 
15 — — — — — — stars T+07 

15 275 10 25 270 15 30 stars T+07 
13 — — — — — — stars FGIO 

10 291 9 10 291 9 10 X-ray HBIO 

5 340 20 10 320 25 15 stars T+07 

19 308 19 19 295 25 25 stars W+09 
10 188 14 14 184 18 18 PNe N+09 
30 — — — — — — stars T+07 

16 491 22 22 538 29 29 X-ray NM09 

20 — — — — — — stars W+09 

10 234 17 19 231 23 26 PNe dL+08 

30 — — — — — — stars T+07 
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Table 1 — Continued 



Name 



Type Mb 



logio 



Vio 



V: 



— V20 + — probe 



ref. 



3203 


so/a 


-19.89 


0.25 


10.47 


229 


7 


7 


3379 


E 


-19.84 


0.11 


10.69 


206 


31 


13 


3957 


so/a 


-19.24 


0.20 


10.38 


199 


13 


13 


4710 


so/a 


-19.10 


0.21 


10.12 


182 


10 


10 


4469 


so/a 


-18.77 


0.17 


10.22 


182 


13 


13 



192 42 19 181 48 



22 



stars W+09 

PNe dL+09 

stars W+09 

stars W+09 

stars W+09 



Note. — The galaxies are sorted by absolute magnitude Mb, which is corrected for Galactic extinction; the quoted errors 
include the statistical uncertainties in distance and photometry. The stellar masses are based on B — V colors, and cor- 
respond to a Chabrier IMF (see text). Circular velocities Vip, V15, and V 2 are measured at 10, 15, and 20 kpc , rcspcc- 
tively, and are in units of km s~^. References: dL+08: Ide Lorenzi et alj ll2008f): dL+09: Ide Lorenzi et al.l ||2009|'): FGIO: 
iForestell fc Geb hardt (20i oh: HBIO: iHuinphrev fc Buot^ 1I2OIOI'); K+OO: iKronawitter et al.( 1I2OOOI'): M+llT iMurphv et al.l 
1I2OIID: MB Ol: MagorriarTfc Ballantvnc ('2001'): N+09: Napolit ano et al.l ||2009D: N+11: iNapolitanoet al.l 1120 11"): NM09: 
iNagino fc M atsushita (2009); OP04: O'SuUivan & Pomrian (2004); O+07: lO'Sullivan et a"Ll l l2007f ); T+07: iThoinaTet al.. (,20071) ; 
W+08: IWeiimans et al.l ||2008|') ; W+09: IWilliams et al.l ||2009D . 
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